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Abstract 

Onset of two- (2D) and three-dimensional (3D) double-component convection due 
to different boundary conditions is studied in a diversely oriented infinite slot with 
broken symmetry between the slot conditions for a component. The main focus is 
on the two compensating background gradients. Different component conditions at 
one slot boundary (the distinction boundary) are considered with such a joint com- 
ponent condition at the other (the similarity boundary) as can be both of the flux 
(v = X = 0) an d of the fixed- value (rj = x = 1) type. Also examined are such compo- 
nent conditions at the second boundary (the inverse boundary) as differ from each 
other inversely to the distinction boundary (r/ = and % = 1). In the horizontal slot 
with inviscid fluid and oscillatory primary instability for r] = x, the most unstable 
wavelength being infinite at r/ = x = is rendered finite by flux-component (solute) 
diffusion at the similarity boundary when rj = x = 1- I n the viscous fluid, however, 
such a diffusion of both components enhances the instability efficiency compared to 
rj = x = 0- With all above types of the broken symmetry, small-amplitude convec- 
tion in viscous fluid remains of an oscillatory nature for any slot orientation other 
than the inversely-stratified horizontal one. For n = and x = Ij such a universality 
also involves various abrupt changes in the marginal-stability curves. These come 
from respectively identified mechanisms of switching between dissimilar oscillatory 
patterns. In inviscid fluid, such changes emerge with the zero instability threshold. 
Some of these abrupt changes give rise to new mechanisms for three-dimensionality 
of the instability. Such a mechanism arises in viscous fluid for rj = x = as well. It 
comes with multiplicity and isolated existence of as well as hysteresis between solu- 
tions of the linear stability equations. Both the hysteresis region and the other above 
abrupt 3D changes are described in terms of analogy between the effect of G = k y /k, 
the ratio of the 2D and 3D wave numbers, and that of a 2D ratio between two grav- 
ity components. Other revealed 3D effects are attributable to new manifestations 
of their general mechanism identified by Tsitverblit [Ann. Phys. 322 (2007) 1727]. 
Apart from dissipation, this mechanism also arises from solute diffusion at the sim- 
ilarity boundary for rj = x = 1 an d from differential gradient diffusion at either 
boundary for n = and % = 1. It also incorporates change of the nature of instabil- 
ity from steady to oscillatory. In the context of the steady linear instability, differen- 
tial gradient diffusion is shown to be more effective at the stress-free slot boundary 
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than at the no-slip one. Although the mechanism of finite-amplitude steady convec- 
tion revealed by Tsitverblit [Phys. Lett. A 329 (2004) 445] is most effective herein 
for rj = x = 1) its manifestation remains well-pronounced for r] = \ < 0.5 as well. 
Relevance of this mechanism to abrupt climate change is thus discussed. 
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1 Introduction 

This work addresses manifestations of a broken symmetry in double-component, 
buoyancy- driven convection resulting from the boundary conditions for one 
component being different from those for the other. Such convection has re- 
cently been identified as a fundamental class of pattern-forming hydrodynamic 
instabilities. The objective of this study is to establish the understanding of 
these instabilities for problems where a major element of previously assumed 
reflection symmetry is absent: the distinction between the components coming 
from one boundary of the fluid domain is not reflected at the other. 

The paradigm of double-component convection in pure fluid first arose in the 
context of conventional double-diffusive convection. This is a class of phenom- 
ena resulting from the effects of unequal diffusion coefficients of two density- 
affecting components |T|2|3] . Among numerous natural science and technol- 
ogy applications of this subject emphasized in its initial reviews [I], particular 
attention has subsequently been focused on small-scale oceanography [5], or- 
dinary evolution of stars [5J, geology [7], geodynamo [SJ, and crystal growth 
[9j. Recently, relevance of double-component convection has also been high- 
lighted for the dynamics of proto-neutron stars during core-collapse supernova 
explosions [ID], as well as for colloidal suspensions [IT], and soap films [12]. 

Since the boundary conditions for one component are generically expected 
to be different from those for the other, the effects of such different com- 
ponent conditions are relevant to all the above areas of application of con- 
ventional double-diffusive convection. However, these effects also apply to 
an eddy-diffusion description of large-scale environmental and turbulent pro- 
cesses, where disparity between the component diffusivities may be practi- 
cally negligible. Such processes range from Langmuir circulations [T5fbl] to 
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the global ocean thermohaline circulation |15|ll6|ll7|ll8|ll9f20f21| and associ- 

ated climate change [22.23.24J. In addition, convective flows are commonly 
used in fundamental studies of transition to turbulence [25] and nonlinear 
pattern formation [26] . 

One major subclass of double-component instabilities arising from the effects 
of different boundary conditions comprises phenomena whose nature is concep- 
tually analogous to the classical double-diffusion [X|2|3] . Differential diffusion 
caused by unequal component gradients forming in perturbed state due to the 
different boundary conditions (differential gradient diffusion) triggers convec- 
tion analogously to the effects of disparate diffusivities. Generalizing the idea 
in [27], such analogy has been introduced in [28 29 3 0|3Tj and scrutinized in 
[32]. 

In particular, the nature of an oscillatory instability highlighted in [27] and 
analyzed in [32] is analogous to that in the diffusive regime of the classi- 
cal double-diffusion [Tf2] . The viscous problem with the stratification inverse 
to that in [27] also gives rise to a mechanism of steady convection [28|29] 
that is conceptually analogous to the finger instability in conventional double- 
diffusive convection [T]. This mechanism generates Langmuir circulations in 
the presence of a stable background density stratification [T3p4] . For the com- 
ponent conditions being different only at one boundary and the other bound- 
ary being infinitely distant, the stratification considered in [281129] has been 
more recently treated in [315] . 

For two horizontal component gradients arising in a laterally heated stably 
stratified slot, the effect of different boundary conditions [30] is also analogous 
to the classical double-diffusion [34] . As in conventional double-diffusive con- 
vection [55], in addition, steady finite-amplitude instability is triggered from 
the state of rest by different sidewall boundary conditions for two compen- 
sating horizontal gradients of the components [31J. Arising without the linear 
steady instability of the conduction state [36] , however, such a finite-amplitude 
manifestation of convection in [31] also exposes an oscillatory linear instability 
whose nature is underlain by differential gradient diffusion [32] . 

Effects of different boundary conditions also extend beyond the instabilities be- 
ing due to differential gradient diffusion. As reported in [37], finite-amplitude 
steady convection arises well before onset of the respective linear instability in 
the viscous version of the problem in [27] . It is then generated by the feedback 
coming from nonlinear Rayleigh — Benard convection, despite the stabilizing 
role of differential gradient diffusion. Potential relevance of such a mechanism 
to abruptly changing global environmental phenomena [15, 16.20,21 22 2 3f24"] 




makes its examination under more realistic conditions particularly important. 
In the above studies of slot double-component convection due to different 
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boundary conditions, the condition for either component at one slot boundary 
has been identical to the respective condition at the other. The considered 
problems have thus been reflectionally symmetric across the slot. Being an 
important initial simplification, the reflection symmetry is however unlikely 
to be maintained in real- world applications of the effects of different boundary 
conditions. This is also particularly relevant if the component-dependent forces 
other than the buoyancy forces are considered, as suggested in [32|37] . 

Manifestation of physical laws in the absence of certain symmetries underly- 
ing them could make both the laws themselves and their broken symmetries 
hardly recognizable. This has been repeatedly illustrated in elementary parti- 
cle physics and cosmological theories of unification of fundamental forces [38] . 
Another illustration is the irreversibility in statistical mechanics [39] and its 
generalized (symmetric) interpretation [10]. In geophysics, the global ocean 
thermohaline circulation also involves asymmetries [TT|18|I19||4"T| . Clarification 
of the nature of such asymmetries is viewed as critical for understanding past 
and predicting future major changes of the Earth climate [22]. 

For classic hydrodynamic instabilities, one can refer to the structure of mul- 
tiple steady flows in the Taylor experiment [12], where translation invariance 
is broken by end walls. Complex as it becomes when the cylinder aspect ratio 
is increased |43j, this structure is not expected to transform into that in the 
translationally invariant problem even when the aspect ratio tends to infinity 
|44j . In addition, if oscillatory instability arises in a flow that is both reflection- 
ally and translationally symmetric, the corresponding Hopf bifurcation would 
give rise to two respectively symmetric oscillatory branches [45] . Elimination 
of one of the symmetries from such a system could thus have a major effect 
on the structure and nature of its nonequilibrium flows. 

The general objective of the present work has been to provide a comprehensive 
insight into the effects of different boundary conditions in a slot where the 
previously assumed boundary conditions symmetry is broken. In a class of such 
problems, the component conditions are different only at one slot boundary 
(hereafter, the distinction boundary). These problems are addressed for such 
a joint component condition at the other (hereafter, the similarity boundary) 
as can range from the flux to the fixed-value type. Being referred to as the 
inverse boundary, this other boundary is also considered with such different 
component conditions as are prescribed oppositely to the distinction boundary. 

Among consequences of the broken symmetry is a universality of oscillatory 
manifestation of the effects of different boundary conditions in viscous fluid. 
For all above types of the boundary conditions, the steady linear instability 
analogous to that in [28|29] transforms into an oscillatory one for any slot devi- 
ation from the horizontal orientation. Three-dimensionality and, eventually, an 
isolated nonlinearity, hysteresis, and other abrupt changes in such oscillatory- 
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linear-instability curves are underlain by the respective horizontal-slot steady 
instabilities. Other new three-dimensional effects come from abruptly emerg- 
ing zero thresholds of the inviscid oscillatory instability. Most pronouncedly 
manifested herein for the fixed-value similarity boundary, the mechanism of 
finite-amplitude steady convection [37] is also relevant when the component 
condition at this boundary is closer to the flux than to the fixed-value type. 



2 The problem formulation and solution procedures 

2.1 The problem and governing equations 

A general case of the considered problem is illustrated in Fig.[U where 9 (> in 
Fig. [T]) is the angle between the direction opposite to the gravity and that of the 
across-slot coordinate axis in an infinite slot with pure fluid. The component 
gradients in Fig.[T]are represented by the Rayleigh numbers Ra = gaATd 3 / 
and Ra s = — gfiips j dx)d A j kv = fiRa. Here, x is the (dimensional) across-slot 
coordinate, d is the width of the slot, AT is the (dimensional) conduction-state 
difference between the values of temperature (the component with the fixed- 
value condition at the distinction boundary) at the boundaries with smaller 
and larger across-slot coordinates, ds/dx is the (dimensional) derivative of 
solute concentration, the component with the flux condition, at the distinction 
boundary, a is the coefficient of thermal expansion, /3 is the coefficient of 
the density variation due to the variation of solute concentration, g is the 
gravitational acceleration, v is the kinematic viscosity, and k = = k>s is 
the diffusivity of both components. The bar means that the respective variable 
is dimensional. Unless explicitly stated otherwise, Ra > and Ra s > as well 
as n = 1 (i.e., the compensating background gradients) are assumed. 

As in [28 29 30 3 T|I32|I37| . the component diffusivities are set equal to eliminate 
the classical double- diffusive effects. Such an approach has also been adopted 
in most studies of conventional double-diffusive convection, where the com- 
ponents with unequal diffusivities were not distinguished from each other in 
terms of their boundary conditions. In principle, equal diffusivities can also 
be experimentally modeled with two solutes [46] . The Prandtl number, which 
would then be significantly different from the present Pr = 6.7, is not expected 
to have a qualitative effect on the main results and physical interpretations 
discussed herein. Equal diffusivities could be interpreted as eddy transport 
coefficients as well, as in [T3~|14|ll9j . Pr = 6.7 is then also within the range of 
realistic values for all diffusion coefficients to be of the eddy type. 

For facilitating comparison of the results for an inclined or vertical slot with 
those for a horizontal slot, this study is focused on the exactly compensating 
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background gradients (/i = 1). This eliminates the along-slot motion arising 
when the slot orientation differs from horizontal. In particular, transformation 
of the oscillatory instabilities at 9 G [0, 7r) into the respective steady instabil- 
ities at 9 = 7r can thus be analyzed in the framework of the effects of differ- 
ent boundary conditions alone. Such compensating gradients have also been 
adopted in many studies of conventional double-diffusive convection [351T4"?] . 

The equations describing the two-dimensional (2D) problem in Fig. [I] can be 
written as follows: 

d( ^ dip d( dip dC, 1 dt ds . ^ 
dr dx dy dy dx Pr dx dx 



1 dt 9s e | d 2 ( | d 2 ( 
Pr dy dy dx 2 dy 



^9+^ + ^, (1) 



= &ip_ + dhp_ ^ 
dx 2 dy 2 ' 

9^ + »±d^_»j L d^ = J_d% + d% . = 12 

dr dx dy dy dx Pr dx 2 dy 2 



Here £i and £ 2 stand for t and s, the across-slot, u, and along-slot, v, velocities 
are 

dip dip 
dy : dx' 

vorticity 

dv du 
dx dy ' 

Pr = u/k is the Prandtl number, r is the time, x G (—1/2,1/2), y G 
(— A/2, A/2), and A = X/d is the specified along-slot period. 

Eqs. P]) — (J3J) are considered along with wall boundary conditions for ( and ip 

C = 7±f^, ^ = (x = ±l/2, -\/2<y< A/2), (4) 



where j± = 1 and 7± = stand for the no-slip and stress-free boundaries, 
respectively, as well as with wall boundary conditions for t and s 

t=™ (o: 1/2, -A/2 <y< A/2), (5) 

— = -//.Ra = -i?a s (x = -1/2, -A/2 < y < A/2, y ^ 0), 
Ra s ds 

(1 - Xo)(s - —) + Xo(^ + ^a s ) = (a; = -1/2, y = 0), (6) 
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dt Rcl 
0--v)(fc +Ra)+ V (t+— )=0 (a: =1/2, -A/2 < y < A/2), (7) 

(9s Ra s 
;i - + #a s ) + X (s + — ) = (x = 1/2, -A/2 < y < A/2, y ^ 0), 

(1 " + ^-) + + = (x = 1/2, y = 0) (8) 



and periodic boundary conditions in the along-slot direction 

{(,,V2)=«.,-A/2), 

(-1/2 < x < 1/2). (9) 

Here £ stands for (, ip, t, and s, Xo — for x < 1 or in all simulations of 
the linear stability for 9 = and Xo — 1 otherwise, xi = 1 f° r 6 = ^/2 and 
Xi = X f° r ^ = 0. In the present study, Pr = 6.7 as well as either rj = x — 
orr/ = x= lor?7 = and x = 1 unless r) = % £ (0, 1) is specified for 6* = 0. 

In (J6j) and (jSJ), specification of the middle values of s identifies the solute scale 
and the solution phase for ?7 = x<lat# = 0, symmetrically fixes the solute 
scale in the temporal simulations fori] = x = Oat 9 = ir/2, and fixes the 
solution phase in such simulations for r/ = x = l a t# = 0. The phase of a 
nontrivial steady solution for 8 = when x = 1 and either rj = 1 or rj = was 
selected by continuation in x and rj from the progenitor of such a solution at 
rj = 1 and x = in [52"f5T] . where the phase and the solute scale were fixed at 
the periodic condition boundaries. The latter approach becomes inconsistent 
with the natural solute scale specification in the nontrivial steady solutions for 
X = 1. It also leads to early spurious oscillatory instabilities of the background 
state for 77 = x = at 6 = 0. It was not thus employed in the reported results. 

In boundary conditions (j3J) — (Q, switching the distinction and the other bound- 
ary designations in Fig. [T] is equivalent to transformation 

(x, y, 9, Ra, Ra s ) T ^ (-x, -y, 9 + n, -Ra, -Ra s ) T . (10) 



Since Eqs. ([T]) — flS} are invariant under QUI) , one can consider only such slot 
orientation to the gravity as is depicted in Fig. [lj 

Discretized by central finite differences, the steady version of Eqs. (JTJ) — 
and boundary conditions (jlj) — © was treated with the Euler — Newton and 
Keller arclength |48j continuation algorithms to trace out bifurcating branches 
[SO] . These algorithms were based on the Harwell MA32 Fortran routine. The 
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along-slot period A = 2 was prescribed. The grid with 33 nodes in the across- 
slot direction was used in the computations [nxx A(nx+1) with nx = 33], as in 
|28f29f30f31f32f37] . As already indicated, temporal behavior of the linearized 
version of Eqs. ([T]) — (J3J) and boundary conditions (j3J — <Q was also examined 
near the onset of oscillatory instability of the conduction state. Such an ex- 
amination was conducted with the implicit method and time step 5r = 0.05. 



2.2 Linear stability calculations 



With the state of rest being the background flow for \i — 1, the Fourier mode 
of a three-dimensional (3D) marginally unstable oscillatory perturbation with 
angular frequency u and a wave number k = [k 2 + k 2 ) 1 / 2 having y and z 
components k y and k z , respectively, can be written as 

[u'(x),t'(x), s >( x )f e i("T±k y y+k z z) + cc _ ( U ) 



(The z-axis is orthogonal to the x-y plane in Fig. [T] and is directed towards 
the reader.) Here [u'(x), f(x), s'(x)] T is the Fourier-mode part depending on 
the across-slot coordinate alone, the prime near a flow variable denotes such 
part in the perturbation of the variable. Expression ffTTl) has been introduced 
into the linearized 3D governing equations that were nondimensionalized con- 
sistently with (2D) Eqs. ([1]) — ([3]) and rewritten in terms of across-slot velocity 
u, temperature t, and solute concentration s. This leads to: 

d 2 d 2 
(^-k 2 )(}- 2 -k 2 - l u J )u = 
dx A dx z 

=F Ra[iGk—{t - s) sin 9 ± k 2 {t - s) cos 6], (12) 
(-^- k2 ~ iu;Pr )& = u>> i = 1 ' 2 ' ( 13 ) 



where G = k y /k, u = —u'Pr, ^ = t — t'/Ra, and £2 = s = s'/Ra s (Ra s = 
Ra). 

The variable u is subject to boundary conditions 



where 7+ = 0, 1 and 7_ = 0,1. Unless the values of r )± are explicitly given, the 
boundaries are assumed to be either both stress-free (7± = 0) or both no-slip 
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(7± — 1)- Such boundary conditions for u are used along with 
ds 

t = ^ = (x= -1/2) (15) 
dx 



and 



;i - ^ + rjb= (1 - x)^ + X s = (x = 1/2), (16) 



where either 77 = x = Oorr] = x = lorr] = and x — 1- 

Let [£t(x),t(x), s(x)] T be a solution of the e i(wr+A:H?/+fcz2) -version of Eqs. ( TT2]) 
and (fT3l) and boundary conditions ( fT4|) — ( fl6l) at some Ra c and u; c for a value 
of # and the orientation of x and y axes as in Fig. [U If (x, y, 6) T is transformed 
into (—x, —y, 6 + ti) t in boundary conditions (fill) — ([IB]) (i.e., if the distinction 
and the other boundary designations in Fig.[T]are switched), [u(x),t(x), s(x)] T 
would then be a solution of the e*^ T ~ fcH?/+fcz2 )-version of Eqs. ffl2l) and (fl3|) at 
—Ra c and the same u c . This symmetry is related to invariance (flOj) of Eqs. 
(II]) — ©• It allows to consider only one orientation of the slot boundaries to 
the gravity, provided that both signs of k y in (II ip are taken into account. 

For examination of the linear stability of the conduction state in inviscid fluid, 

d 2 d 
iu(— - k 2 )u = Ra[±iGk—(t - S) sm9 + k 2 (t- s) cosfl] (17) 
dx 2 dx 

was used along with 
d 2 

— -k 2 -iu)^ = u, z = l,2, (18) 



1 dx 



where the Rayleigh numbers are defined as Ra = gaATd 3 /n 2 and Ra s = 
—gP(ds/dx)d 4: /K 2 = Ra. Here u is nondimensionalized with n/d 2 as opposed 
to vjd 2 in Eqs. (Tl2|) and (JT3J) . The definitions of Ra, Ra s , and u for inviscid 
fluid are then different from the respective definitions for viscous fluid. Such 
definitions are thus used below according to the type of fluid in question. 

The boundary conditions for inviscid fluid are 

u = ( x = ±l/2) (19) 

along with ( Fl5l) and ( TT6l) . 
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For 7] = and x — 1; Eqs. ( T12"]) and (TT3]) along with boundary conditions 
( IT^|) — ( TT6|) for 7_ = 7+ as well as Eqs. ( ITT)) and (fl8|) and boundary conditions 
f )15p . f)16p . and ( I19p are invariant under transformation 

[u(x),t(x), s(x), Of ^ [u(-x), S(-X),t(-X), 7T - fl] T . (20) 

For i] = and x = 1, therefore, 9 G [0, 7r/2] and 9 G [vr/2, 7r] are equivalent in 
the context of (1201) . For Ra ^ Ra s at 9 = and 9 = tt as well as for 7_ 7^ 7+ 
in viscous fluid, such an equivalence also implies additional transformations. 
These are specified when the respective results are discussed below. 

The inviscid problem in a horizontal slot was also examined for the Ra and 
Ra s being independent of each other. In particular, 

d 2 

-k 2 )u = k 2 (Rat- Ra s s) cos 9 (21) 

ax 2 

was considered along with Eqs. (fl8l) and boundary conditions ( fT5l) . ( [161) . and 
f fl9]) . This was done at given values of Ra s and 9 = for any of the pairs of 
7] and x specified above. As already indicated, such a problem for rj = and 
X = 1 is equivalent to that with the fixed values of Ra at 9 = 71. 

For fixed k y and i?a c and u c are found by searching in the Ra-u domain for 
the smallest Ra at which the complex matrix resulting from the application of 
boundary conditions f JT4|) — f fT6|) to the general solution of Eqs. fTl2|) and ffT3]) is 
singular. The same procedure is applied to the general solution either of Eqs. 
f JTTj) and f fTHj) or of Eqs. ffl8]) and (1211) for 9 = with boundary conditions 
( fl5l) . ( fl6|) . and ( fl9l) . NAG Fortran routines were employed for this purpose. 

Once Ra c {ky) and uj c {ky) have been found for a given /c° and fixed k z , the 
corresponding values of these parameters, Ra c {k y ) and u c (k y ), at a nearby 
fcy are computed by the Euler — Newton continuation method. The latter is 
applied to the solution of equation 

F[Ra(k y ),u(k y ),k y ,k z } = 0, (22) 

where F[Ra(k y ),u(k y ), k y , k z ] stands for the (complex) determinant of the 
matrix resulting from the application of boundary conditions ( TT^j) — f JT6|) or 
(JTSJ), ( ITS"]) , and f TTU]) to the general solution of the respective set of differen- 
tial equations. Due to the use of standard Fortran routines, the Jacobian of 
{Re[F(Ra,u, k y , k z )}, Im[F(Ra,u, k y , k z )]} T (with respect to Ra and uj) and 
dF(Ra,u, k y , k z )/dk y were computed with numerical differentiation. When 
3D effects were anticipated, they were studied by repeating the procedure just 
described for different k z > 0. 
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Since the linear instability to steady disturbances in viscous fluid was found to 
arise only in a horizontal slot, the corresponding problem is discussed in the 2D 
framework alone. For the compensating background gradients, in particular, 
co = and k y = k are set in Eqs. (TT2~|) and (1T31 to obtain 

d 2 d 
(—-k 2 ) 2 u = -Ra\ik—(t - s) sin9 + k 2 (t - s) cos9], (23) 
dx l dx 

^-k% = u, i = l,2. (24) 



At 9 = tt, the marginal-stability boundaries described by Eqs. f l23|) and (1241 
and boundary conditions (|14p — ffTBl) were found for any of the above pairs of 
i] and x- F° r V — and \ — 1 3 they were found at 6* = as well, in view of 
(1201) . For 6 = and 9 = it, steady linear instability was also examined for the 
Ra and Ra s being independent of each other: 

d 2 

(— - k 2 ) 2 u = -k 2 {Rat- Ra s §) cos 6 (25) 



and Eqs. ([241 were used either at # = ir for any of the pairs of r] and \ specified 
above or at 9 = for 77 = and \ = 1 alone. 

The general solution either of Eqs. ([231 and ([24"]) or of Eqs. (|24|) and (1231) 
is obtained analytically. Boundary conditions ( 1T4|) — (1T61) are then applied to 
such a general solution. In the former case, the smallest Ra, Ra c (k), at which 
the resulting matrix becomes singular is searched for at different k (and found 
only for 9 = tt, besides 9 = when 77 = and x — !)• m the latter case, such 
a search yields either the smallest Ra s , Ra s c (k), for a fixed Ra when 9 = 7T or 
the smallest i?a, Ra c (k), for a fixed i?a s when # = at 77 = and \ = 1- 



3 Inviscid fluid 



5J ^ = 



5 JJ 77 = X = 

When = 0, the anticipated linear instability for 77 = x is oscillatory [2~Tf32] . 
The corresponding inviscid linear stability problem is generally described by 
Eqs. (fl~8|) and ([2T]) and boundary conditions (1T5|) . (1T61) . and (1T9|) . For flux com- 
ponent conditions at the similarity boundary, the marginal-stability curves, 
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Ra c (k) and oj c (k), are illustrated in Fig. [2(a). The qualitative similarities be- 
tween them and such curves in the symmetric (77 = 1, x = 0) case [Fig- 2(a) 
in [32]] are the result of the same basic physics of the instability. This physics 
and its implications in Fig. [2(a) are summarized just below. As explained in 
|32j . one can consider only a standing- wave perturbation [45J. 

In the end of a rotation cycle of a perturbation cell, a potential energy of com- 
ponent perturbation stratifications is generated. Due to differential gradient 
diffusion at the distinction boundary, it is utilized by the cell in the beginning 
of the cycle of rotation in the opposite sense. The maximal amount of such 
energy depends on the instability horizontal scale. This scale determines the 
time available for the two components of a fluid element to diffuse. Ra c (k) 
thus decreases with the increase of the horizontal wavelength and becomes 
minimal for a given Ra s as k — > [Fig. [2(a)]. 

For efficient utilization of the potential energy, the frequency with which the 
marginally unstable cells change their sense of rotation also has to resonantly 
match the time for vertical diffusion naturally specified by the instability wave- 
length. As a consequence, u c (k) — > as k — » and u c (k) grows with the 
increase of k from [Fig. [2(a)]. 

As the wave number is further increased, however, the wavelength time for 
the manifestation of differential diffusion eventually becomes insufficient for 
the cell oscillation amplitude to grow. To afford more time for such diffusion, 
u) c (k) thus decreases when certain values of k are exceeded [Fig. [2(a)]. This 
leads to an inconsistency between the diffusion time afforded by the oscillation 
frequency and that specified by the instability wavelength. Due to efficiency of 
the instability mechanism being then reduced, the instability fails to develop 
above a critical value of the wave number. With the increase of Ra s and 
resulting growth of Ra c , the enhanced gradient disparity in the perturbed 
state intensifies the energy transfer to the perturbation cells. The minimal 
unstable wavelength thus decreases as Ra s grows [Fig. [2(a)]. 

Compared to 77 = 1 and x = 0, however, differential gradient diffusion arises 
at 77 = x = only from one boundary, the distinction boundary. More time is 
thus required for such process to be effective as that at 77 = 1 and x — 0- For 
this reason, u c (k) are smaller at 77 = x = [Fig- 12(a)] than at 77 = 1 and % = 
[Fig. 2(a) in [32] ] for all k > 0. The restriction of differential gradient diffusion 
for i] = x — to a vicinity of the single boundary also explains why the 
respective intervals of unstable k terminate at slightly smaller wave numbers 
compared to rj = 1 and % = 0. For the same reason, Ra c (k) are higher at 
77 = x — than at rj = 1 and x = f° r a U k outside a vicinity of k = 0. 

When the instability horizontal scale is large enough, however, differential 
gradient diffusion becomes as effective for development of the instability at 
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one boundary as at both. Additional temperature diffusion due to the fixed- 
value condition at the second boundary for 77 = 1 and x = then becomes a 
stabilizing factor compared to 77 = x — 0. It reduces the ability of unstable 
temperature stratification to generate the potential energy of solute perturba- 
tion stratification. In the vicinity of zero wave number, therefore, Ra c (k) are 
smaller for 77 = x = than for 77 = 1 and x = 0. For the Ra s in the present 
Fig. |2f a) and in Fig. 2(a) of [32J, this takes place when k < 0.2 for Ra s = 5000, 
k < 0.1 for Ra s = 10000 and Ra s = 20000, and k < 0.06 for Ra s = 50000. 

As for rj = 1 and x = 0, the fact that the most unstable wave number for 
77 = x = is zero makes the determination of exact values of Ra c (0) and 
group velocity w£(0) = du c (0)/dk relevant. Using the same long- wavelength 
expansion as employed for 77 = 1 and x = in [32J, one obtains 

<4(0) = ^® = u = v/i^/12 (26) 
and then 

Ra c (0) = (2Ra s + 5040)/156. (27) 

These expressions differ from the respective expressions for 77 = 1 and x = 
[32] only by the denominator in Eq. (127]) . The numerical data underlying the 
marginal-stability curves in the present Fig. [2(a) accurately coincide with Eqs. 
(j2BD and fl27]) . For the two compensating gradients, Ra c (0) = Ra s , 

Ra c (0) = 2520/77, tu c k (0) = co = ^210/77. (28) 
5.i.jg r/ = x=l 

For fixed- value component conditions at the similarity boundary, the marginal- 
stability curves are illustrated in Fig. 12(b). They are qualitatively different 
from such curves both for the symmetric case [Fig. 2(a) in |32J] and for the flux 
similarity boundary [Fig. [2(a) herein] by the relative stability of the vicinity of 
zero wave number. For any Ra s , Ra c (k) increases abruptly when k decreases 
below a certain value, for which Ra c (k) is minimal. An immediate vicinity of 
k = is also stable for any Ra. 

The abrupt increase of Ra c (k) with decreasing k is associated with solute 
diffusion at the similarity boundary. Such diffusion diminishes the role played 
by the potential energy of solute perturbation stratification in the instability 
mechanism. When the instability wavelength is relatively short, however, the 
enhancement of differential gradient diffusion with growth of the wavelength 
still lowers the Ra c (k). 



13 



Below a certain k, however, neutralization of the solute perturbation scale 
forming at the small Ra c (k) by its diffusion at the similarity boundary becomes 
more important than the effect of differential gradient diffusion. Ra c (k) then 
grows with the wavelength, to generate a higher solute perturbation amplitude 
at the same Ra s . Arising from a fixed Ra s , however, the solute perturbation is 
merely eliminated by its diffusion at the similarity boundary for any Ra when 
the wave number reaches an immediate vicinity of k = 0. The instability of 
such small k thus fails to develop [Fig. EJ^b)]. 

3.1.3 7] = and x = 1 

For Eqs. (ITS]) and (12Tj) . the present results would also be applicable to 9 = tt 
when transformation fl20]) is accompanied by Ra c > Ra s c and Ra s t— > Ra. To 
avoid confusion, however, they are discussed only in terms of 9 = 0. 

With the inversely different component conditions at the boundaries, a hori- 
zontal slot combines elements of both steady and oscillatory instability mech- 
anisms. For 9 = 0, in particular, the component stratifications near the dis- 
tinction boundary give rise to a mechanism of oscillatory instability of the 
type discussed above. 

Near the inverse boundary for 9 = 0, however, the component stratifications 
correspond to a mechanism of steady instability of the type reported in [2%|T25] . 
Also arising from differential gradient diffusion, the mechanism of such an 
instability leads to amplitude growth of only such perturbation cells as do 
not change their sense of rotation. This mechanism is expected to affect the 
oscillatory instability coming from the distinction boundary. 

The perturbation cells whose sense of rotation changes periodically in time 
would however be located closer to the distinction boundary. Generated by 
such cells, the disparity between component perturbation gradients near this 
boundary is expected to exceed the opposite one near the inverse boundary. 
For small and intermediate wavelengths of the oscillatory perturbation, there- 
fore, the time available for differential gradient diffusion would prevent such 
a (relatively small-gradient-disparity) process near the inverse boundary from 
being effective. Largely specified by the distinction boundary, Ra c {k) in Fig. 
Ws) thus decrease with k for such wavelengths. 

Above a critical wavelength, however, the time available for differential diffu- 
sion becomes sufficient for such a process near the inverse boundary to notice- 
ably damp the oscillatory perturbation. The potential energy of component 
perturbation stratification generated by differential gradient diffusion at the 
distinction boundary is appreciably reduced by such a process at the inverse 
boundary. For this reason, Ra c (k) increases with k decreasing below a certain 
value [Fig. 12(c)]. 



14 



When the wavelength becomes long enough, the combined effects of differ- 
ential gradient diffusion near the distinction and inverse boundaries result in 
larger Ra c (k) corresponding to the smaller Ra s [Fig. |2fc)]. Ra s is not only a 
measure of stability near the distinction boundary. Via the inverse boundary, 
it also controls the steady opposition to growth of the oscillatory perturba- 
tion. In particular, this opposition is the stronger the more time is afforded 
for differential gradient diffusion with the growing wavelength. As the relative 
disparity between the diffusion times [represented by the respective l/u c (k), 
Fig.|2Jc)] for different Ra s also increases with the wavelength, such effect of the 
inverse boundary becomes more important than the effect of the distinction 
boundary. Larger Ra s are thus destabilized by the smaller Ra c (k). 

In the immediate vicinity of zero wave number, the (large) time available for 
differential gradient diffusion makes the effect of such a process near the inverse 
boundary comparable with that near the distinction boundary. This takes 
place despite the smaller gradient disparity formed near the former boundary. 
As a consequence, the oscillatory instability fails to develop [Fig. 12(c)]. 

3.2 9e[0,n/2] 

3.2.1 General 

Oscillatory linear instability in a vertical slot with viscous fluid is discussed 
in Sec. l4.1.2.2l below. The broken boundary conditions symmetry then results 
in only one traveling wave. Such a disturbance arises from differential gradi- 
ent diffusion, whether due to the distinction boundary alone or to both the 
distinction and the inverse boundaries. The direction of propagation of such a 
traveling wave matches the sign of the background contribution to the density 
(Fig. [1]) at a sidewall with different boundary conditions from a (component) 
variable whose flux is prescribed there by dHj) — (|HJ). [This also applies both to 
the single traveling wave for rj = and X = 1 and to either of the symmetri- 
cally counter-propagating waves for rj = 1 and x = (Fig. 7 of [32])] F° r the 
configuration in Fig. [T] (6 = 7r/2), such a traveling wave has +k y in (jllj) . 

Growth of the disturbance with — k y in ( II ip is inconsistent with the effect of 
differential gradient diffusion combined with the along-slot gravity component. 
2D instability to the standing-wave disturbance at 9 = cannot thus transform 
into the respective instability to the e l ^ T ~ kvy ' '-mode traveling wave when 6 is 
increased from 0. As a consequence, instability to the latter traveling wave 
was found to vanish precipitously when 9 increases from for all types of the 
boundary conditions. 

The vanishing 2D instability is then replaced by the respective instability to 
3D perturbations, since the latter perturbations are less sensitive to the above 
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asymmetry introduced by the along-slot gravity component. (In particular, 
the along-slot gravity does not affect a perturbation with k y = and k z ^ 0.) 
Largely driven by the across-slot gravity, such a 3D instability also has to 
vanish as 9 increases further. In view of these considerations, it is only the 
e i(ur+k y y+k z z)_ moc \ ie instability that is discussed herein below both for inviscid 
and for viscous fluid. 

For any type of the boundary conditions, the inviscid e^ bJT+kyy+k:iZ ' > -mode in- 
stability with certain k z > possesses intervals of k y for which Ra c (k y ) = 
Lo c {k y ) = at some 9 G (0, 7r/2], as discussed below. Such a zero u c , however, 
does not necessarily mean that the instability is steady. In particular, steady 
instabilities with finite Ra c are absent in viscous fluid for 9 G (0, 7r) and the 
boundary conditions considered herein. The zero-o; c instability arising say at 
9 = 7r/2 from the inviscid equations for oscillatory marginally unstable per- 
turbation then has no immediate connection with the instability of a steady 
origin. It seems thus reasonable to expect that the perturbation developing at 
any small Ra(k y ) > for some 9 G (0, tt/2] and k z > would generally have a 
small u)(k y ) > as well. This is what is implied below, although such an effect 
of along-slot gravity is referred to as direct. 



3.2.2 r] = x = 

At 9 = 7r/2, the (2D) e l< ^ r+fc!/y )-mode instability would generally have to be 
manifested in the form of traveling cells whose sense of rotation is constant 
in time. In inviscid fluid, such a manifestation also has to take place at any 
Ra > 0. In the absence of dissipation, even an infinitesimal horizontal density 
difference resulting from differential gradient diffusion could drive the pertur- 
bation. Eq. fTl7|) then implies that such instability with Ra c (k y ) = also has 
to have ui c {k y ) = 0. For 9 G (0,7r/2), however, the effect of across-slot gravity 
opposes amplitude growth of the perturbation cells whose sense of rotation 
does not change (Sec. 13. 1.1)) . 

The combination of the effects of across-slot and along-slot gravity thus results 
in a decrease of u c {k y ) with 9 increasing from to 7r/2 [Fig. [3]^a) and (b)]. 
This allows part of the rotation energy of a perturbation cell to come from 
the energy directly contributed by the along-slot gravity component in the 
current cycle of rotation. In the end of a cell rotation cycle, the whole rotation 
energy is transformed into the potential energy of perturbation stratification 
due to the across-slot gravity component. This potential energy is released in 
the next rotation cycle. As 9 grows from 0, the Ra c (k y ) thus also decrease. 

At the longest 2D wavelengths for 7] = x — 0, the across-slot gravity compo- 
nent acts most effectively in opposing the direct contribution of the along-slot 
component to the rotation-energy increase. On the other hand, the shortest 



16 



scales are least effective in accommodating the contribution of the along-slot 
gravity. It is thus a set of intermediate 2D wavelengths that becomes most 
unstable when 9 = n/2 is approached [Fig. Efb)]. Decreasing for all k y as 
9 — > 7r/2, such Ra c (k y ) and cu c {k y ) also become identically zero at 9 = ir/2. 

Upon introduction of the long-wavelength expansion used above [32] into Eqs. 
ffTTl) and ([18]) with k = k y , the k\ order of the e i(ajr+fcl/?/) -version of these 
equations and boundary conditions f [T5l) . f[T6]) . and ([19]) (r] — x — 0) yields 

fta c (0) = 24^/(2 cos 9 + u sin 9), (29) 

which is consistent with Eq. ( 1261) for 6> = 0. The values of Ra c (0) and cjfc(0) = 
Uq estimated from the numerical data underlying the (2D) marginal-stability 
curves in Fig. [3](a) and (b) [as well as from such data for k y < (for which 
w < as well) not reported herein] were found to be fairly consistent with 
Eq. PD- 

3.2.3 v = X = l 

At # = 7r/2 for 77 = % = 1, the inviscid e^^^^-mode instability also has 
to arise at any Ra(k y ) > 0, and thus u c (k y ) = for such instability as well. 
When only the across-slot gravity component is present (9 = 0), the vicinity 
of zero wave number is most stable to the disturbances that change their sense 
of rotation periodically in time [Fig. Ht^a)]. It is therefore this region of k y that 
least opposes the direct contribution of along-slot gravity to the cell rotation 
energy. As 9 increases from 0, 2D Ra c (k y ) thus decreases most significantly 
near k y = [Fig.H^a) — (e)]. Closer to 9 = ir/2, however, the Ra c {k y ) decreases 
noticeably and tends to as 9 — > n/2 for all k y [Fig. HJf)]. 

As discussed in [32J, three-dimensionality of most unstable disturbances could 
be a consequence of two general mathematical conditions. One of them (con- 
dition I) is dependence of Ra c only on wave number modulus k, as in 

Ra c (k y , k z ) = Ra c [(k 2 y + k 2 z ) l '\ 0], (30) 

at a single value of some parameter (9 = in this case) in whose vicinity Ra c 
depends on both components of k. The other condition (condition II) is the 
existence of an interval where Ra c (k) is growing with decreasing k at this value 
of the parameter. Three-dimensionality of the instability in a vicinity of the 
above parameter value then follows from the assumption that these conditions 
would largely apply at nearby values of such a parameter as well. 

As discussed in Sec. 13.1.21 above, condition II holds at 9 = [Fig. BJa)] near 
k = due to solute diffusion at the similarity boundary. This process is 
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therefore responsible for the three-dimensionality of instability in Fig. H^b) 
and (c). With the Ra c (k y )\ 9=7T / 2 = 0, the fast decrease of 2D Ra c (k y ) near 
k y = with 9 growing from prevents such 3D disturbances from being 
dominant for the larger 9. That no 3D most unstable disturbances are found 
in Fig. [3]^c) is thus a consequence both of condition II being then not met for 
9 = [Fig. [3](a)] and of the absence of another mechanism for 3D instability. 

3.2.4 %D disturbances for 77 = and x = 1 

3.2.4.1 General. Transformation (1201) makes the results of this Sec. 13. 2.41 
also applicable to 9 6 [tv/2,tv]. Assuming ( 1201) . they are however discussed 
below only in terms of 9 G [0, 7r/2]. As seen from Fig. |5l any small increase of 
9 from leads to abrupt changes in the 2D marginal-stability curve for 77 = 
and x = 1. These changes are associated with emergence of an interval of k y 
for which Ra c (k y ) = u c (k y ) = 0. Such interval is born with its width and the 
lower limit tending to zero as 9 — > 0. Both these parameters of the interval 
then become finite when 9 increases from 0. They also continue to grow with 
9 increasing further (see the solid lines in Fig. [6]). 

In viscous fluid, steady instability for 9 = arises at finite viscous Ra c (k) 
(Fig. Cj). (It precedes the respective instability to standing wave.) The inviscid 
steady Ra c {k) are then zero for all k > [Eq. fl23l) for 9 = 0]. As discussed 
below (Sec. 14. 1.31) . however, such 2D steady viscous instability transforms into 
an oscillatory instability to traveling wave when 9 is increased from 0. As 9 > 0, 
therefore, the argument just used does not apply. For 9 > 0, the inviscid 2D 
instability with Ra c {k y ) = uj c {k y ) = still ought to arise from the type of 
perturbation that is dominant at 9 = tt/2, where the largest zero-threshold 
interval forms [Fig. HD^f)]. Generally, this has to be a traveling wave (whose 
speed turns infinitesimal with u c = 0). Such a perturbation is also the first to 
become unstable in viscous fluid for 9 e (0,7r/2]. Its manifestation would be 
most convenient to analyze in the presence of along-slot gravity alone. 



3.2.4.2 Effect of the along-slot gravity. In the presence of along-slot 
gravity alone [Fig. 0(f) for k z = 0], Ra c (k y ) = u c (k y ) = for k y > ky l (0) fa 
2.4677 (Table[T]for 9 = 7r/2 and k z = 0). Just below this value of k y , however, 
both Ra c {k y ) and u c (k y ) abruptly increase to finite magnitudes. Additional 
calculations for k y > 6 and 9 = n/2 also suggest that the upper limit of the 
interval with zero Ra c {k y ) and u c {k y ) is at infinity. 

As discussed below for viscous fluid (Sec. 14. 1.2.21) . differential gradient diffu- 
sion for f] = and x — 1 ac ts in concert at the vertical-slot boundaries. It gives 
rise to horizontal density differences that favor growth of the e^'^^-mode 
traveling wave. In particular, the density perturbation generated around a 
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clockwise-(counterclockwise-)rotating cell is largely specified by negative (pos- 
itive) solute perturbation at the distinction boundary and temperature pertur- 
bation at the inverse boundary. This intensifies such downwards-propagating 
small- amplitude convective cells with a steady sense of rotation. 

In the absence of dissipation, the horizontal density differences just described 
would give rise to the instability even when they are infinitesimal. Such an 
inviscid vertical-slot instability can thus arise for any Ra > no matter how 
much the effect of differential gradient diffusion diminishes with decreasing 
the wavelength. This is seen in Fig. Eff) (k z = 0) for k y > fcjj^O) « 2.4677. 

For 7] = and x = 1, however, either component forms a diffusion gradi- 
ent at one of the boundaries. Such a diffusion is also the more effective the 
longer the wavelength is. When the wavelength exceeds a critical magnitude, 
therefore, the perturbation of either component could be neutralized by its 
diffusion at one of the boundaries. The infinitesimal horizontal density dif- 
ferences arising just above Ra c (k y ) = 0J c {k y ) = are thus eliminated by the 
component diffusion. Both Ra c (k y ) and u c (k y ) must then increase from zero 
as k y decreases below the critical value [k y = k y l (0) ~ 2.4677]. This increase is 
precipitous since such growing u c {k y ) also lowers the efficiency of differential 
gradient diffusion, due to inconsistency between the u c (k y ) and the increasing 
wavelength. 



3.2.4.3 Combined effects of the along-slot and across-slot gravity. 

The oscillatory instability for 9 = arises only at finite Ra c (k y ) (Fig. [5]). 
It is characterized by a standing wave, i.e. by convective cells whose sense 
of rotation changes periodically in time. Development of such a perturbation 
is inconsistent with that of a traveling wave arising due to the along-slot 
gravity, for such traveling-wave convective cells do not change their sense of 
rotation. For 9 G (0, 7r/2), the latter perturbation is also favored and opposed 
by the effects of across-slot gravity at the inverse and distinction boundaries, 
respectively, particularly when it is manifested with Ra c (k y ) = u c (k y ) = 0. 
The steadily rotating cells would thus have to be localized near the inverse 
boundary to the extent 9 is close to 0. 

Due to the expansion of the (steadily rotating) convective cells towards the 
distinction boundary with 9 growing from to n/2, however, an opposition 
at this boundary to the steady sense of cell rotation is generally relevant for 
any 9 < n/2. Its relative role in rotation of a convective cell depends on the 
instability wavelength. As such wavelength decreases, in particular, the effect 
of across-slot gravity at the distinction boundary (opposing the steady sense 
of cell rotation) becomes relatively more pronounced with respect to that at 
the inverse boundary (favoring the steadily rotating cells). In addition, the 
shorter the wavelength the smaller the relative portion of streamline particles 
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with across-slot density differences compared to that with such along-slot dif- 
ferences. The decreasing (increasing) wavelength thus also enhances (reduces) 
the effect of across-slot gravity with respect to that of along-slot gravity. 

The minimal wavelength above which the zero-threshold steadily rotating cells 
are dominant then depends on 9. Increasing from with 9 decreasing from 7r/2 
[Fig. [6] for k z — 0], it tends to infinity as 9 — > (Fig. [5]). Just above such a 
critical k y [= k 2u (0) (Table Q]for k z = 0)] for 9 G (0, vr/2), the steadily rotating 
perturbation cells fail to grow at infinitesimal Ra. The Ra c {k y ) and u c (k y ) 
thus increase from 0. 

As such Ra c {k y ) and u c {k y ) grow from 0, however, the effect of across-slot 
gravity becomes increasingly more relevant. The efficiency of differential gra- 
dient diffusion for the steadily rotating perturbation cells then decreases. The 
Ra c (k y ) and co c (k y ) thus grow precipitously to such k y . Relatively localized 
near the inverse boundary (to the extent 9 is close to 0), the steadily rotat- 
ing convective cells then also transform into such perturbation cells as are 
relatively localized near the distinction boundary and change their sense of 
rotation with an adequate frequency. Control over the instability disturbances 
is thus largely transferred to the across-slot gravity component, whose action 
still remains affected by the direct contribution from the along-slot component. 

When 9 = n/2 and k y < k y l (0)\e =7T /2 ~ 2.4677, an infinitesimally small Ra{k y ) 
is neutralized by diffusion and thus fails to generate the instability. This takes 
place when such a diffusion process is equally active at both boundaries, as 
is expected at 9 = n/2. The closer 9 is to 0, however, the more localized the 
steadily rotating perturbation cells are near the inverse boundary. Compared 
to 9 = 7r/2, the overall effect of diffusion for such cells when 9 G (0,7r/2) 
thus becomes asymmetrically divided between the components. Temperature 
diffusion at the distinction boundary is then less active than solute diffusion at 
the inverse boundary. For this reason, ky l (0) [below which the Ra c (k y ) abruptly 
increases from 0] also decreases with 9 G (0,7r/2), from ky l (0)\ e=7r / 2 ~ 2.4677 
to as 9 -> (Fig. |5]and Fig. |6]for k z = 0). 



3.2.5 3D disturbances for r] = and x = 1 

3.2.5.1 General. Although transformation ( 120]) makes the results of this 
Sec. 13.2.51 also applicable to 9 G [ir/2,ir], except for Sec. I3.2.5.5[ they are 
discussed below only in terms of 9 G [0,7r/2]. For any 9 G (0,7r/2) in Tabled] 
and Fig. [6J the 3D lower and upper limits of the above zero-threshold interval 
of k y , k 2 y \k z ) and k^(k z ), decrease from kf° = fcj(0) and k 2 y u0 = fc^(0), 
respectively, with k z increasing from 0. They also continue to decrease with 
k z increasing further. For 9 = tt/2, this formally applies only to ky l (k z ), since 
ky U (k z ) is then at infinity for any finite k z > 0. In addition, another zero- 
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threshold area arises from the vicinity of k y = when k z > [see the dashed 
lines in Fig. Efa) — (e)]. However, the lower and upper limits of the latter area, 
ky l (k z ) and ky U (k z ), respectively, increase with growing k z (Tabled! 9 ^ tt/2). 

At a certain k z > 0, the interval of finite Ra c (k y ) and 0J c (k y ) between ky U (k z ) 
and k 2l (k z ) thus vanishes, due to these parameters of k y merging with each 
other. A single continuous interval of zero Ra c (k y ) and u c (k y ) is then formed 
[the dotted lines in Fig. H^a) — (e)]. Eventually, such a continuous interval 
also vanishes when ky l (k x ) and ky U (k z ) merge at a still larger k z , leaving only 
nonzero Ra c and uj c [the dash-dot lines in Fig. [Bfa) — (e)]. Quantitative details 
of the behavior just described are reported in Table [H 

When 3D perturbations in Fig. [6] are dominant, their behavior cannot be 
explained only in the framework of the scenario emphasized in Sec. l3.2.3l above. 
Other mechanisms causing such a three-dimensionality would thus also have 
to exist. These mechanisms are associated with the nature of the 2D interval 
with Ra c (k y ) = u) c {k y ) = discussed in Sec. 13.2.41 above. 



3.2.5.2 k 2 y \k z ) and k 2u (k z ). For a given 9 G (0,vr/2) and k z = 0, the 
2D Ra c (k y ) and oo c (k y ) increase from zero either when k y decreases below ky 10 
or when it increases above k 2u0 . This takes place because for the respective 
wavelength intervals (X y > 2n/k 210 and X y < 2n/k 2u0 ) at the fixed 9 e (0, vr/2), 
the perturbation cells with a steady sense of rotation cannot be destabilized 
by the effects of differential gradient diffusion at infinitesimal Ra (see Sec. 
13.2.41 above). Depending on the relative roles of the along-slot and across-slot 
gravity components, ky 10 and k 2u0 are also the closer to the closer 9 is to 0. 

Independent of the orientation of the axis of rotation of a convective cell, the 
effect of across-slot gravity does not change when k z increases from 0. In this 
context, one could therefore refer to ky 10 and k 2u0 as the respective critical 
values of k. With k z growing from 0, however, the (relevant) projection of 
along-slot gravity on the axis orthogonal to the axis of cell rotation decreases. 

Indeed, Eqs. (|T7|) and (fTB|) for k = (k y , k z ) are identical to these equations 
for k = (k y ,0) with that sin# in the latter is replaced by (k y /k)sm9 in the 
former. The relative role of along-slot gravity is thus diminished with respect 
to that in the 2D problem for k = (k y ,0). This is similar to decreasing 9 
below the considered value. With increasing k z , therefore, the actual values 
of ky\k z ) and k 2 y u (k, z ) must be smaller than q y l {k z ) = [{ky 10 ) 2 — k 2 ] 1 / 2 and 
q 2u (k z ) = [(k 2u0 ) 2 - A^] 1/2 , respectively, as in Tabled] for any 9 ^ vr/2. 

For 9 = it/2, the relative role of the along-slot gravity component with respect 
to that of the (absent) across-slot component would remain infinite for any 
finite k z . This has to lead to k 2l {k z ) = q 2l {k z ), as is seen from Table[T](# = tt/2). 
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In addition, the infinite value of ky U ° is retained both by qy U {k z )^ k >Q and by 
k 2 y U (h) lkz>Q [Fig. Mf) and Table ffl(0 = tt/2)]. 

For 6* = 7r/2, ky l {k z ) and ky U {k z ) should thus be exactly specified by fc^ zo and 
fcy u0 (= oo). These latter have to serve as the critical values of k. Indeed, Eqs. 
(fTTD and (HHP for = vr/2 suggest that 

Ra c (ky,k z )k y /k = Ra c (k,0). (31) 
This means that Ra c (k y , fc z )i e=7r / 2 = -<=>• Ra c (k, O)| 0=7r ^ 2 = for 

*r lte , /2 <(^+« ,/2 <*r lfc , /2 =~. (32) 

which is a condition for k[— (ky + fc^) 1 ^ 2 ] alone. 

For 9 = 7r/2, 3D disturbances are thus most unstable at least for k y be- 
tween qf!(k z ), n ;o (< k^°, a .„) and A;,^ ,. ,„. This could cause such a three- 
dimensionality for small tt/2 — 9 as well. For any 9 G (0,7r/2), however, 
ky l {k z ) < qy l {k z ) (because G = k y /k < 1), as discussed above. The latter 
inequality is also an independent cause for dominance of 3D disturbances 
at 9 G (0,7r/2). Additional causes of perturbation three-dimensionality for 
9 G (0,7r/2) are associated with another zero-threshold area in 
plane. This area is described in Table [TJby ky l (k z ) and ky U (k z ). 



3.2.5.3 General on ky l (k z ) and ky U {k z ). The nature of the area with 
Ra c (k y ,k z ) = u c (ky,k z ) = whose boundaries are described in Table [TJ 
by ky l (k z ) and ky U (k z ) is fundamentally three-dimensional. First note that 
Eq. ([30]) is satisfied by Ra c (k y , k z )\ d€{0j7T/2 ) for k y = 0: Ra c (0, k z )^ Q ^ /2) = 
i?a c (0, fc)|0 6 ( O = Ra c (k, 0)| e=0 / cos 9. In a vicinity of k = [Fig. [Sj, there- 
fore, Ra c (0, fc) |g g( - / 2 ) decreases with increasing Conditions I and II dis- 
cussed in Sec. 13.2.31 are then met. A region of small k z is thus expected to be 
dominated by 3D most unstable disturbances with small G = k y /k > 0. This 
observation, however, does not unravel the behavior of ky l (k z ) and ky U (k z ). 

Eqs. dUD and ([18]) for k y > also suggest that kf(k z ) . e=e ^ and k^ u (k z ). d=dg 

specify such limits of the interval of k over which Ra c (k, 0) = u c (k, 0) = for 
9 = 9^ as are equal to the respective limits of two 2D intervals of k y over which 
Ra c (k y , 0) = u c (k y , 0) = for two smaller values of 9. These are such respective 
9 = 9 2 {< 9 3 ) as tan6» 2 = [h y (k z )/k 3 (k z )] tan0 3 , where k 3y (k z ) = kf(k z )^ =0 ^ 

k y U ( k *)\e=e 3 and MM = [ k U k ^ + ^ ■ For ^ fixed 9 E (M/2), &e 
increase of G = k y /k and k z from could thus give rise to a 3D area with 
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Ra c (k y , k z ) = (jJ c (k y , k z ) = analogously to the emergence of the 2D interval 
between k 210 and k 2u0 near k y = when 9 increases from (Fig. [5]). 

Three-dimensionality of the perturbation allows the relative effects of the two 
gravity components to be varied independently of 9, and in particular to mim- 
ick the vicinity of 9 = at any 9 G (0, tt/2). As parameters of k y between which 
Ra c (k y , k z ) = u c (k y , k z ) = 0, k^{k z ) and ky U (k z ) thus not only depend on the 
relation between the gravity components. They also specify this relation. 



3.2.5.4 k^(k z ) and k* u (k z ) for k z 0. For a fixed 9 G (0,vr/2), the 
effects of across-slot and along-slot gravity are represented in Figs. |5] and EJ^f), 
respectively. In view of Eq. (1301) . the dependence on k y for k z = in Fig.[5](# = 
0) also represents the dependence on k for k z > 0. The dependence on k z > in 
Fig. ©(f) is then given by Eq. When k -> 0, Re[du(k, 0)/dx] x u c (k, 0) x 
1 both at 9 = and at 9 = tt/2 [as implied by the scaled continuity and Figs. 
and 0(f)] . From Eqs. ffTTl) and ffT8l) . the general solution for u(x) then suggests 
that Ra c (k, 0)\ 6=0 x (l/k 2 ) and i?a c (/c, 0) | 0=7r / 2 x (l/k) as A; — > 0. Defining 

= iMMV^/iMM),^, therefore, ££ /2 (&)| fe ^ Q x fc(-> 0). 

For (k z /k y ) < oo [k z = 0(k y )] when A; — > 0, therefore, Eq. (I3TT) implies that the 
effect of along-slot gravity dominates that of across-slot gravity. For a finite 
9 G (0,7r/2), this would be inconsistent with the existence of an infinitesimal 
area in the k y -k z plane (for infinitesimal k) where Ra c (k y , k z ) = u c (k y , k z ) = 0. 
Such an area could arise only when e = {k y /k z ) — > [i.e., k y = o{k z )} as well. 

When e along with k and i? c (&) = Rl /2 (k)/k, the necessary balance 
between the effects of the two gravity components [analogous to their 2D 
balance given by < (k y / tan6>) < oo as k y (= k 2u0 ) — > and 9 — > 0], 

< lim [{k 2 /k y )R n /2 {k)/k} = lim[R c {k){k 2 /ky)} < oo, (33) 

e->0 c->0 
k^0 fc->0 



implies that < ik 2 jk y ) < oo for k — > and e — >■ 0. In particular, numerical 
computations show that R c (k) is maximal at k — >■ and that 

lim RJk) » 1.055. (34) 



Using the l'Hospital rule, one could then obtain, in particular, 

(k 2 /k y ) ~ (fc'/ig ~ fc/e ~ ~ 2*/(^) ~ 2*,/(^) 
dk z d(^k )-i^/9 r^(^2:)ii/2 dky dJvy 



Vi- 



de de 'die 2 )' ^(e 2 ) 1 die 2 ) ll dik 2 ) 
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where ~ denotes the asymptotic equivalence. With such necessary conditions, 
the relation between the effects of the two gravity components is similar to 
that in the 2D problem when 9 — > 0. For any finite 9 G (0, vr/2), such infinites- 
imal ky l (k z )[= o(k z )} and ky U {k z )[= o{k z )} as maintain f l35j) could thus arise, 
Ra c (k y , k t ) = u c (k y , k t ) = for k y (k z ) G [tff(k t ), k% u (k z )}. 



3.2.5.5 ky l (k z ) and k y u (k z ) for finite k z . Equivalence relations ( )35|) im- 
ply that for k y = k 3y (k z )[= k y n {k z ),ky U {k z )\, dk y /dk z — > as k z — » and 
e — > 0. k 3y (k z ) then has a local minimum in this limit, since d 2 k y /dk 2 z > as 
k z — > and £ — > for such a k y {k z ). One could thus expect that the respective 
dk y /dk z be positive within a small (positive) vicinity of k z = 0. It is generally 
shown below that when k y = k 3y (k z ), dk y /dk z > for any k z (> 0). 

As discussed in Sec. 13.2.41 ky and k 2u0 grow with increasing (decreasing) 
9 G (0, 7r/2) [# G (7r/2, 7r)], when the along-slot gravity is enhanced with 
respect to the across-slot gravity. Eq. (TIT|) also suggests that the variation of 
G(k, k x ) = {k 2 — k 2 z ) x l 2 jk causes qualitatively the same effect on the relation 
between the gravity components as that of | tan#|. k 3 (k z ){= \k 2 y {k z ) + k 2 ] 1 / 2 } 
thus has to grow with G: dG(k, k z )/dk > for k = k 3 (k z ). This yields 

— < — — > — (36) 
dk 3 k 3 dk z k z 



Since e — > as k z — > for k y — k 3y (k z ), one could expect that G(k 3 ,k z ) = 
kz y /k 3 [~ k y /k z = e for infinitesimal k z and k y = k 3y {k z ) = o(k z )) would 
increase with k z within a small vicinity of k z = 0. Such a behavior is found 
to hold for any finite k z (> 0) as well: (13"B"1) also implies that dG(k, k z )/dk z > 
for k = k 3 (k z ). With G = k y /(k 2 + k 2 z ) x l 2 and k y = k 3y (k z ), however, 



This means that dk 3y {k z ) / dk z > k 3y (k z )/k z > 0. 

ky l (k z ) and ky U (k z ) thus have to only increase with growing k z . Such a behavior 
is consistent with the data in Tabled] Due to k 2l {k z ) and k 2u {k z ) decreasing 
with growing k z , it also explains both the merging of ky U (k z ) with k 2l (k z ) and 
that of k^(k z ) with k 2u (k z ). 
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Another aspect of the data in Table[T]is that for a fixed k z , k y l (k z ) and k y u (k z ) 
are the closer to the closer 9 is to tt/2. This is associated with their ratio be- 
tween the gravity components, their G(k y , k z ) tan 9, being a function of k alone, 
as suggested by Eq. (IT7j) . Such a function thus has to be proportional to 1/k. 
For k z (> 0) being fixed and k y = k 3y (k z , 9), therefore, d[k y (k z , 9) tan 9] /d9 = 



implies dk y (k z ,9)/d9 < [dk y (k Zl 9)/d9 > 0] for 9 e (0,tt/2) [9 e (tt/2, tt)]: 



Since ky l (k z ) and k^ u (k z ) increase with 9 approaching n/2, the ultimate van- 
ishing of all areas with Ra c (k y , k z ) = u c (k y ,k z ) = 0, when k y l {k z ) and k y u (k z ) 
merge, takes place at the larger k z the closer 9 is to tt/2 (Tabled]). For 9 = tt/2, 
the absence of across-slot gravity makes the relation between the gravity com- 
ponents independent of G. Eq. fl3Tl) thus prohibits the existence of ky l (k z ) ^ 
and k y u {k z ) ^ for any k z . 



4 Viscous fluid 

4-1 Small- amplitude oscillatory convection 
4.1.1 9 = 

4.1.1.1 Diffusion for f] = x- Despite dissipation, different component 
conditions at one boundary give rise to oscillatory instability at 9 = both 
for 7] = x — and for r\ = \ — 1- For the no-slip boundaries, this is illustrated 
in Fig. [HI As could be expected, Ra c (k) and u c (k) are higher for both values 
of r] = x [Fig. [8](a)] than for r\ = 1 and x = [Fig. 5(a) in [32]]. Fig. Mjo) also 
exhibits decaying oscillations of SRa c (k) = Ra c (k)^ =x=0 — Ra c (k)^ v=x=1 and 
Suj c (k) = u c (k) lri=x=0 - uj c (k) {v=x=1 . 

The behavior of 5Ra c (k) and 5co c (k) in Fig. [S](b) arises from a more effective 
role of the distinction boundary for 77 = \ — 1- Compared to 77 = x = 0, 
diffusion at the similarity boundary reduces the component perturbation scales 
for 77 = x = 1 by the same fraction (that increases with the wavelength). 
Specified by the background scales, the gradient disparity at the distinction 
boundary then grows with respect to such a reduced component scale. The 
relative amplitude of convective motion thus also grows, generating a relatively 
larger gradient disparity for the next rotation cycle. With such more efficient 
mechanism, the small-A; instability for n = x — 1 precedes that for 77 = x = 0- 

Figs. [3(a), (f) and [T0l(a).(f) illustrate the stage of potential energy release via 



dk y (k z ,9) 
09 



sin(2fl) 



2 



ky(k z ,9). 



(38) 
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differential gradient diffusion. The ratio of the streamfunction perturbation 
scale to that of either component for 77 = x — 1 [Fig- flOT a) and (f)] exceeds 
this ratio for 77 = x — [Fig. [H^a) and (f)]. Thus relatively more intensive at 
77 = x — I; such a convective motion gives rise to new component perturba- 
tion stratifications [Figs. Mjo) — (d) and [101(b) — (d)]. The latter arise from the 
respective background gradients with the ones for 77 = % — 1 being smaller. 

Via the higher gradient of temperature diffusion at the distinction boundary, 
the perturbation stratification thus formed opposes the current sense of cell 
rotation [Figs. [9](c),(d) and fT0l(c).(d)]. Due to the higher efficiency of such 
differential diffusion, this opposition is more pronounced for the fixed-value 
similarity boundary than for the flux one. With respect to the perturbation 
scale of either component, in particular, the streamfunction scale for 77 = x = 1 
[Fig. [roT c).(d)] becomes smaller than that for r\ = x = [Fig. |^c),(d)]. The 
relatively unequal potential energies of perturbation stratification so generated 
are then utilized by the reversely rotating cells [Figs. Ete),(f) and [TOf e).(f)]. 



4.1.1.2 Diffusion and dissipation for 77 = x- A consequence of differ- 
ential gradient diffusion being more efficient for 77 = x = 1 than for 77 = x — 
is thus a relatively larger variation of the respective velocity scale. For about 
the same variation time at a fixed wavelength, this implies steeper spatial 
velocity gradients for 77 = x = 1 and thus the respectively greater dissipation. 

Of two dissipation mechanisms affecting the instability [32], one is merely asso- 
ciated with damping all motions. Its overall effect is the greater the shorter the 
instability wavelength is. The other mechanism causes an efficiency reduction 
for the instability feedback. Such feedback links the component perturbation 
potential energy generated when a cell rotates in one sense and the intensity 
of rotation of such a cell in the opposite sense. The effect of this mechanism 
depends on the along-slot part of dissipation of a convective cell. It is the 
increasing role of the latter dissipation mechanism that causes Ra c (k) to rise 
infinitely with k decreasing to for both rj = x- The higher efficiency of 
differential gradient diffusion for 77 = % — 1 thus becomes relevant for such k. 

When the wavelength increases, the growing destabilizing contribution of dif- 
fusion at the similarity boundary is increasingly opposed by the respective 
enhancement of only the second of the two above dissipation mechanisms. 
Both these counter effects are commensurately augmented by the growing 
wavelength. The enhancement of overall dissipation for rj — % — 1 (with re- 
spect to rj = x — 0) would thus have to remain of a limited relative significance 
compared to the respectively higher efficiency of differential gradient diffusion. 

For sufficiently long waves, therefore, it is the efficiency of differential gradient 
diffusion that specifies at which value of 77 = x the instability sets in first: 
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5Ra c (k) > [Fig. E^b)]. For a fixed wavelength, in addition, the optimal 
frequency with which the convective cells change their sense of rotation has 
to be largely specified by the background gradient of the stably stratified 
component. For /i = 1, therefore, 5u c (k) > for such k. 

When the instability wavelength decreases, the effect of the first dissipation 
mechanism is enhanced, due to the across-slot motion being augmented. The 
role of the overall dissipation disparity between r\ = x — 1 an d r\ = X — 
then grows with respect to the (diminishing) effect of diffusion at the similarity 
boundary. \SRa c (k)\ and |5a; c (/c)| thus decay with increasing k [Fig. Efb)]. 

However, the dissipation enhancement for rj = x — 1 eventually dominates 
the respectively higher efficiency of differential diffusion underlying it. 5Ra c (k) 
thus becomes negative [Fig. E^b)], and then so does Su c (k). The overall effi- 
ciency of the combined effects of diffusion and dissipation then becomes higher 
for rj = x = 0. An additional dissipation arising above the k where SRa c (k) = 
would therefore be greater for rj = x — than for r\ — % — 1. When the rel- 
ative role of the disparity between such additional dissipations is augmented 
sufficiently with k increasing further, the sign of 5Ra c (k) changes again. 

The oscillatory nature of decay of SRa c (k) and 5u c (k) in Fig. MJo) is thus 
associated with the additional dissipation arising above a critical value of k 
where 5Ra c (k) = acting against the increase of \SRa c (k)\. The value of 
7] = x & t which the overall efficiency of diffusion and dissipation becomes 
higher above such a critical k is also characterized by more dissipation. 

The values of k where 5u) c (k) = in Fig. 0(b) slightly exceed the respective k 
where 5Ra c (k) = 0. A comparatively more time for diffusion and dissipation is 
thus provided for the value of 7] = x that has just [when 5Ra c (k) = 0] become 
a less efficient combination of these processes. Over a short interval of k, this 
outweighs the mismatch between the signs of 5u c (k) and 5Ra c (k). 

Unreported herein, the values of k at which 5Ra c (k) and 5u c (k) change sign 
for stress-free boundary conditions were found to be respectively smaller than 
those for the no-slip conditions. This could be associated with the higher no- 
slip Ra c (k) being more important for delaying the effects of dissipation than 
the decrease of (along-slot) dissipation due to the stress-free boundaries. 



4.1.1.3 r] = and \ — 1. Oscillatory instability for r] = and x = 1 a t 

9 = {6 = 7r) arises where the fluid is already unstable to steady disturbances 
(Fig. [7j). It is however instructive to comment on its Ra c (k) and oj c {k) com- 
pared to those for i] = x- For rj = and x = 1> the oscillatory perturbation 
has to be localized near the distinction (inverse) boundary. With respect to 
t] = x, this implies a higher efficiency of differential gradient diffusion at this 
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boundary. The relative portion of moving fluid particles with the generated 
horizontal density differences is larger for 77 = and x = 1 than for 77 = x- 

For i] = and x = 1> however, differential gradient diffusion at the inverse 
(distinction) boundary opposes amplitude growth of a convective cell whose 
sense of rotation changes periodically This opposition is dominant for rela- 
tively long wavelengths, whereas the higher efficiency of oscillatory across-slot 
motions is more important for the short waves. In viscous fluid, such a higher 
efficiency at the large k for r\ = and x = 1 is also expected to generate a 
relatively higher dissipation. This would however be offset by the across-slot 
cell path being shorter for rj = and x = 1 than for 1] = x- The effects of 
dissipation thus have to be relatively little important in the present context. 

As k exceeds a certain value, Ra c (k) for rj = and x = 1 thus changes from 
being larger to being smaller than that at either of 77 = x- F° r inviscid fluid 
and for viscous fluid with no-slip and stress-free boundary conditions, such a 
value is between 4.9 and 5.3. For u c (k), such a value is between 5.6 and 5.8. 
[The data for most such cases are depicted in Figs. [S^a) and@Ja), partly in 
Fig. [5j as well as in Figs. [7]and[8]^a).] The larger critical values of k for ui c {k) 
further relatively diminish the effect of differential gradient diffusion at the 
inverse (distinction) boundary. Until then, this effect remains more important 
than the match between u c (k) and the respective stable solute stratification. 



4.1.2 9e[0,n/2] 

4.1.2.1 General. The marginal- stability curves in Fig. [8]^a) meet condi- 
tions I and II emphasized in Sec. 13.2.31 above. Here condition II holds due to 
the second dissipation mechanism (Sec. 14. 1 . 1 .2j) for both rj = x and also due 
to solute diffusion at the similarity boundary for rj = x = 1- 3D disturbances 
are thus the first to arise near k y = for small 9 > [Figs. [TlTa) and [T27 a)]. 
Their dominance intervals of k y are shorter than those for 77 = 1 and x = 
[32] . This results from a less favorable combination of the effects of two gravity 
components. The 2D instability for rj = 1 and x = se ^ s in at 6 1 = before 
and at 9 = it /2 after those for 77 = In particular, the pattern of two counter- 
traveling waves for rj = 1 and x = 0at# = 7r/2 (Fig. 7 in [32J) generates more 
dissipation than such a single-wave pattern for 77 = % considered below. 

As for 77 = 1 and x = 0, however, the intervals of k y with 3D most unstable 
disturbances vanish for 77 = x only when 9 = 71 /2 (Figs. [TT1 and [T2l) . With 
Re[du(k)/dx]\0 =Q x u c (k)\ e=0 x 1 and Re[du(k,0)/dx]\g =7T / 2 x u c (k, 0)| 9=7r/2 x 
1 as k — )► [Figs. Eta) and fTBT a)]. the general solution for u(x) from Eqs. f|T2l 
and ([I3D suggests that (Sec. 13.2.5.41) [R^ /2 (k)/G]\ k ^ x (k/G)\ k ^ > 
G\t^o — 0(k)(—t 0). Arising from the across-slot gravity alone, the three- 
dimensionality is thus retained by a sufficiently small G = k y /k so long as 



28 



9 < 7r/2. It could also be viewed as coming from conditions I and II at G = 0. 
As k y — > for a fixed k z > 0, in particular, the Ra c (k y , k z ) cos 9 and co c (k y , k z ) 
tend to the respective 2D Ra c {k z ) and co c (k z ) at 9 = [Figs. ISlfa). ITT] and[T2"]. 

In view of transformation (1201) . Ra c (k y , k z ) and u c (k y , k z ) for 77 = and x = 1 
when 9 G [0,7r/2) are addressed when G (7r/2,7r] is considered in Sec. 14.1.31 
below. Since Eq. ( 1HT|) also applies to viscous fluid, the 2D disturbances must 
be most unstable at 9 = tt/2 for any 77 and This is seen from Figs. ITBT a) and 
IT4T a). The former figure is independent of 77 = x- The respective instability 
mechanisms for 77 = x thus have to be such as are essentially due only to the 
distinction boundary. By virtue of Eq. ( )3~T]) . this could be discussed in the 2D 
framework alone. The mechanisms of 2D instability at 9 = 7r/2 for 77 = x are 
thus considered below along with such a mechanism for 77 = and X = 1. 



4.1.2.2 2D instability mechanisms for 9 = tt/2. The patterns illus- 
trated in Figs. [TSl [161 an d EH are traveling waves propagating in the gravity 
direction. Their nature can be clarified if one assumes that in the marginally 
unstable state, the speeds of their propagation adequately match the respec- 
tive convective velocities oriented downwards. In the reference frame moving 
with such a pattern, the component perturbations are then transported ver- 
tically only upwards. This is the direction where the convective velocities are 
augmented by the moving reference frame. 

The slot area behind such a propagating clockwise-(counterclockwise-)rotating 
cell is supplied with the temperature and solute perturbations from the region 
near the left (right) sidewalk The background component values are relatively 
high (low) there. Occupied by a counterclockwise-(clockwise-)rotating cell, 
such an area is thus largely characterized by the positive (negative) component 
perturbations (Figs. [151 (El and [T7|) . 

Due to differential gradient diffusion, the density near the distinction boundary 
to the left is specified mainly by the solute perturbation there. At the similarity 
boundary to the right, however, the temperature and solute isolines behave 
identically to each other for either 77 = x (Figs. [TBI and [TBI) . The density excess 
there thus has to be zero. (Due to the fixed solute sidewall values at y — 0, this 
holds only approximately in Fig. [13) For 77 = and x — 1 (Fig. fTT]) . differential 
gradient diffusion at the inverse boundary to the right results in the density 
there being specified mainly by the respective temperature perturbation. 

A propagating convective cell that rotates counterclockwise (clockwise) is thus 
characterized by positive (negative) horizontal density differences between the 
streamline regions at the left and right sidewalls. This is what drives convective 
motion for such a cell. For 77 = x (Figs- EH an d [TBI) , such differences arise 
mainly from the positive (negative) solute perturbation at the left sidewall and 
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zero density perturbation at the right boundary. For rj = and x = 1 (Fig- 
[17)1 . they are mainly due to the positive (negative) solute and temperature 
perturbations at the left and right sidewalls, respectively 

Such a downwards-propagating counterclockwise-(clockwise-)rotating connec- 
tive cell also has its component perturbations near the left (right) sidewall 
practically steady. The convective-cell velocities there are largely offset by the 
speed of propagation. There is thus a steady horizontal density difference that 
maintains the downwards propagation of such a flow pattern. For 77 = x (Figs. 
[15] and [16]) , this difference is specified by the positive left-sidewall solute per- 
turbation and zero right-sidewall density perturbation for a counterclockwise- 
rotating cell. For 77 = and x = 1 (Fig. [17]) . the steady density difference comes 
from the cells rotating in both senses. It is due to the positive (negative) left- 
(right-) and zero right- (left-) sidewall solute (temperature) perturbation for a 
counterclockwise- (clockwise-)rotating cell. 



4.1.2.3 Marginal-stability curves. For any pair of r] and x just con_ 
sidered at 9 = ir/2, the propagating 2D convective pattern gives rise to such 
a distribution of the component perturbations as favors its convective motion 
and maintains its direction of propagation. This is implemented due to differ- 
ential gradient diffusion. Such a process takes place at the distinction sidewall 
alone for f] = x an d at both vertical boundaries for 77 = and x — 1- 

For any rj and Xi l ne instability mechanism at 9 = tt/2 is underlain by hor- 
izontal density differences accompanying the flow pattern. The independence 
of 2D Ra c {k y ) and ui c {k y ) in Fig. [TBTa) of the value of i] = x is therefore just 
a manifestation of the component disparity at the distinction boundary being 
unaffected by the orientation of component isolines at the similarity sidewall. 
Eq. f]3T]) then suggests that such independence is retained for k z > as well. 
Since the increase of k z for a fixed k y decreases the overall wavelength, u c (k y ) 
in Figs. [TBT a) and [T4T a) also increases with k z . 

Convective cells driven by such across-slot density differences between their 
vertically moving particles make their along-slot dissipation a part of the insta- 
bility feedback. Analogously to the second dissipation mechanism for standing- 
wave perturbation at 9 = 0, discussed in Sec. 14.1.1.21 above, this causes an 
infinite growth of 2D Ra c (k y ) with k y decreasing to at # = 7r/2. At = 0, 
such a growth for rj = x = 1 is also caused by solute diffusion at the similarity 
boundary. For rj = x, the 2D marginal-stability curves at 9 = [Fig. E^a)] are 
thus smoothly transformed into those at 9 = tt/2 [Figs. [TTIT2] and [TBT a)]. For 
7] = and x = 1 at = tt/2, the increase of 2D Ra c (k y ) near k y = is also due 
to neutralization of either component by diffusion. Such a marginal-stability 
curve for 9 decreasing from tt/2 to is discussed in Sec. 14. 1.3.31 below. 
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Inversely affecting the component perturbations at the vertical sidewalls, dif- 
ferential gradient diffusion for 77 = and x = 1 still leads to the respective 
horizontal density differences augmenting each other by their superposition. 
The instability for these 77 and x [Fig- El( a )] thus sets in substantially before 
that for 77 = x [Fig- [131(a)] . Such combination of the effects of differential diffu- 
sion at the boundaries also results in the convection pattern (Fig. [T7|) without 
a slope between its across-slot motion and the horizontal axis. 

Such horizontality of across-slot motion, however, makes the overall cell path 
shorter for 77 = and x = 1 than for 77 = x- With respect to the relative 
streamfunction amplitudes, therefore, the maximal horizontal density differ- 
ences in Fig. [T7] exceed those in Figs. [15] and [161 The intensity of convection 
in the marginally unstable state has to be matched by the speed of pattern 
propagation and thus by the associated u c (k y ). This explains why u c {k y ) are 
respectively smaller for 77 = and x = 1 [Fig- HU(a)], where convective motion 
is relatively weaker, than for rj = x [Fig. [TBT a)]. Also consistent with this 
argument is that the difference between such u c (k y ) grows with k y . 

4.1.3 9e[ir/2,ii} 

4.1.3.1 General on 2D disturbances. For 9 = it, the across-slot grav- 
ity is manifested as a steady instability [28.29]. For the values of 7/ and x 
considered in this study, such an instability is discussed in Sec. 14.2.21 below. 
Also characterized by steadily rotating convective cells, the 2D traveling waves 
arising at 9 = n/2 thus have to transform into the respective steady distur- 
bances when 9 increases above 7r/2. The speed of flow-pattern propagation, 
u c (k y )/k y , then has to decrease to for all k y . For the reflectionally symmetric 
pattern of two counter-propagating waves arising at 9 = tt /2 when 77 = 1 and 
X = [32], such a decrease begins with certain k y at 9 ~ 1.57r/2. 

For the present values of 77 and x-, however, the mechanisms described in Sec. 
14.1.2.21 above have to retain a flow-pattern propagation for any 9 G [ir/2,ir). 
The decrease of u c (k y ) to for any k y (> 0) can thus take place only when 
9 = it, as in Fig. [T37 b) and (c) and in Fig. [T47 b). The phases of nonzero 
component perturbations in Figs. HU [161 an d EH also nearly coincide with the 
streamfunction phases of the opposite sign. This is due to the traveling nature 
of such convective patterns (Sec. 14.1. 2.2p . When u c (k y ) turns with 9 reaching 
tt, however, these relative phases have to become shifted by a quarter of the 
wavelength with respect to each other, as in Fig. 2 of [29] and in Fig. [TBTa). 



4.1.3.2 2D disturbances for rj = x- For 9 = tt/2, the along-slot part of 
dissipation reduces efficiency of the instability feedback. This does not apply 
to 9 = tt, where only the first of the two dissipation mechanisms accentuated 
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in Sec. 14.1.1.21 above is relevant. In the framework of this first mechanism, 
however, the across-slot part of dissipation for 77 = x is a l so larger at 9 = tt/2 
than at 9 = tt. This is due to the slope of such an across-slot cell motion at 9 = 
tt/2 (Figs. [TS] and US}. For r] = x = 0, 2D Ra c (k y ) is thus higher at 9 G [tt/2, tt) 
than at 9 = 71 for any k y [Fig. ITBT a) and (b)]. The difference between such 
Ra c (k y ) also has to be infinite at any 9 G [tt/2, tt) as k y [— o(tt — 9)]— > for any 
uj c (k y )^ k ^ in Fig. [TBT b). due to the effect of along-slot dissipation. Compared 
to u c (k y )^^ — > 0, however, such effect is moderated when u c (k y )^^ x 1. 

As 9 exceeds tt/2, in particular, two joined marginal-stability branches with 
U} c(ky)\k -^0 ~~ ** ^ (due to the across-slot gravity) also isolatedly arise from 
k y = [Fig. ITBTb). 9 = \Att j2\. Their higher Ra c (k y ) has the slightly higher 
u c (ky). Both their Ra c (k y ) are smaller and larger than that of the main branch 
[whose u c (k y )\ k ^ x 1, due to the along-slot gravity] at some k y > and at 
k y — > 0, respectively. [Ra c {k y )^ k ^ Q — > 00 is assumed for the smallest u c (k y ) as 
well.] Growing with 9, the higher- and lower-w c branches with uj c {ky)^ k ^ —> 
meet the main branch at a finite k y and unfold with its smaller- and larger-fc y 
intervals, respectively [Fig. [TBTb). 9 > 1.5n/2]. The branch with u c (k y )^^ Q x 
1 then exists only below such a finite k y as decreases to with 7r — 9. 

For T) = x — 1, solute diffusion at the similarity boundary increases 2D Ra c (k y ) 
infinitely as k y decreases to at 9 = n [Fig. [TBTc)] (Sec. 14.2.2.21 below). At 
9 = tt/2, such Ra c {k y ) is still independent of 77 = % [Fig. ITBTa)]. Its infinite 
increase for k y — > is then due only to the above role of along-slot dissipa- 
tion. When k y decreases, however, the relative portion of streamline particles 
whose horizontal density differences drive a convective cell grows at 9 = 71 / 2 
compared to 9 = 71. This efficiency factor eventually dominates the discrep- 
ancy between the growing stabilizing effect of along-slot dissipation and that 
of solute diffusion at the similarity boundary. The 2D Ra c (k y ) near k y = 
thus becomes smaller at 9 = n/2 [Fig. d3J(a)] than at 9 = n [Fig. E^c)]. With 
9 increasing from tt/2, therefore, the Ra c (k y ) grows for such very small k y and 
decreases elsewhere to transform into its values at 9 = tt. 

4.1.3.3 2D disturbances for rj = and % — 1. For rj = and x — 1) the 

across-slot gravity when 9 = tt (9 = 0) favors and opposes growth of steadily 
rotating convective cells at the distinction (inverse) and inverse (distinction) 
boundaries, respectively. Its effect is manifested in terms of respective along- 
slot density differences arising between the streamline particles moving across 
the slot close to these boundaries. For 9 G (tt/2, tt) [9 G (0, vr/2)], it is combined 
with the effect of along-slot gravity. The latter effect is manifested (Fig. [T71) via 
across-slot density differences between the streamline particles moving along 
the boundaries. As 2D k y — > 0, then, Ra c (k y )\Q =7T /2 < Ra c (k y )\g =n fl [Fig. [T4T b)]. 

Since the instability for 9 = tt/2 is also due to differential gradient diffusion at 
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both sidewalls, 2D Ra c (k y ) for such 9 is smaller than that of steady instability 
for 9 = 7r [9 = 0) so long as the wavelength is large enough for the effect of 
diffusion to be dominant [Fig. Hlf b)]. At 9 = n (9 = 0), however, the steadily 
rotating convective cells are localized near the distinction (inverse) boundary 
[Fig. [T8T a)]. The across-slot part of their dissipation thus decreases compared 
to 9 = 7r/2. With increasing k y , the relative role of the streamline particles 
with along-slot density differences is also enhanced with respect to that with 
across-slot density differences. For sufficiently large k y (k y > 5.3), therefore, 
the Ra c (k y ) for 9 = n (9 = 0) is smaller than that for 9 = tt/2 [Fig. Efb)]. 

As 9 changes from tt/2 to tt (0), the lost contribution of along-slot gravity to 
the intensity of the steadily rotating cells is replaced by the mutually opposing 
effects of across-slot gravity at the distinction and inverse boundaries. This 
weakens such convective motion. To maintain the remaining effect of along- 
slot gravity, the speed of propagation of the flow pattern has to match the 
relative intensity of convection. u c (k y ) then decreases. This shifts the relative 
phases of component and flow perturbations with respect to each other. The 
portion of streamline particles with favorable across-slot density differences 
thus decreases, and the relative effect of along-slot gravity weakens further. 
The role of this feedback depends on the wavelength and 9. 

If the wavelength is too short for retaining a sufficiently long slot interval 
where the across-slot density differences matter, an efficient utilization of the 
along-slot gravity becomes impossible. As a consequence, the perturbation 
of such a wavelength occupying the whole width of the slot cannot persist. 
It then gives way to a convective pattern driven mainly by the across-slot 
gravity. Localized near the distinction (inverse) boundary, such a pattern has 
differently behaving Ra c {k y ). These just exceed the Ra c {k y ) at 9 = tt {9 = 0) 
for the relatively large k y in Fig. ITW b). Mainly underlain by the (steady) effect 
of across-slot gravity, such a pattern is also characterized by ui c {k y ) [the large- 
ly ujciky) in Fig. [IW b)] that are substantially smaller than those at 9 = tt/2. 

For sufficiently large \9 — tt/2\, comparatively abrupt changes in dRa c /dk y and 
dbj c jdk y are distinguishable in Fig. [TU^b). Such changes are a manifestation of 
the transition from the (relatively long-wavelength) convective pattern largely 
driven by the along-slot gravity to the localized (relatively short-wavelength) 
pattern mainly driven by the across-slot gravity. The longer the wavelength is 
(the larger is the ratio between the portions of streamline particles with across- 
slot and along-slot density differences) the more capable its convective pattern 
is of accommodating the effects of across-slot gravity without destroying the 
mechanism by means of which the along-slot gravity drives such convection. 
The abrupt changes thus arise at the smaller k y the more \9 — tt/2\ exceeds 0. 
Their k y in Fig. [TJJb) also tends to zero with 9 approaching tt (0). 
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4.1.3.4 3D disturbances for r\ = x = 1« For 77 = % = 1, solute diffusion 
at the similarity boundary at 9 = tt results in the Ra c (k) for steady insta- 
bility rising to infinity with k decreasing to [Fig. ITBT c)]. as also discussed 
in Sec. !4.2.2.2l below. Conditions I and II (Sec. I3.2.3P for three-dimensionality 
of the instability at small tt — 9 > are thus met. That the instability at 
9 < tt is oscillatory is consistent with the corresponding mechanism for its 
three-dimensionality. Such oscillatory instability with small u c > is a per- 
turbation of the steady instability (u c = 0) at 9 = tt that preserves the effects 
of conditions I and II. 

3D oscillatory disturbances (k z > 0) are thus most unstable near k y = for 
small tt-9 > [Fig.H^a)]. A sufficiently small G = k y /k in Eqs. (pD and flB]) 
also makes the 3D effect of across-slot gravity dominant in the respectively 
small vicinity of k y = so long as 9 7^ tt/2 [Fig. [T9f b).(c)]. Indeed, with 
Re[du(k,0)/dx]\ e=n „ o ( k ) x u c (k, 0), g o(k) x 1 and Re[du(k, 0)/dx]\ e=7T / 2 x 
u c (k,0) le=7r/2 x 1 as k -> [Fig. [TS^a) and (c)], G\ k ^ = 0(k)(-^ 0) 

[Rn-o(k)(k) / G}\k^o x [k/G]\ k ^o > 0. The three-dimensionality for any 9 G 
(tt/2, 7r) could be viewed as coming from conditions I and II at G = as well. 
As k y — > for a fixed k z > 0, in particular, the Ra c (k y , k z ) \ cos#| tends to the 
2D Ra c (k z ) at 9 = n [Figs. HS](c) and US] with the u) c (k v ,k g ), k -)■ (Fig. 
[T9l) . This effect vanishes only when = n/2, where Eq. fl3T|) applies. 

4.1.3.5 3D disturbances for = x = 0. For rj = x — 0, a vicinity of 
fc y = in Fig. [2D] is still dominated by 3D disturbances, despite condition II 
being not met at 9 = 71 [Fig. [TBT b)]. The role of 3D disturbances with small k y 
in Fig. [20] is also enhanced when 9 decreases from tt. Growing with tt — 9, the 
effect of along-slot dissipation on the feedback efficiency increasingly heightens 
the 2D Ra c (k y ) in the vicinity of k y = 0. The three-dimensionality then comes 
from a much more unstable behavior of such Ra c (k y ) at 9 = tt, since the latter 
behavior could be mimicked by Ra c (k y , k z )\g € ( n / 2j - K ) for small G = k y /k. 

For any k z > 0, there have to exist such k y as make G so small that the 
along-slot gravity in Eqs. (fl2l and ([TBI be negligible. The axes of cell rotation 
are then nearly parallel to the along-slot gravity. For these k y , u c (k y ,k z ) and 
)\cos9\ at 9 G (7r/2, 7r) are respectively approximated by and 
2D Ra c {k)\ e=v . In Fig. [TSTb). 2D Ra c (k)^ e=7T reaches its minimum, 720, as 
k p2M9] (Sec. 14.2.2.11 below). For such G and small enough k z , the 
Ra c (k y , k z ) ~ i?a c (0, k z ) ~ Ra c (k z , 0) ~ 720/ 1 cos^. It is thus smaller than the 
smallest 2D Ra c (k y ) {k ^ (-> 00) at G (7r/2,7r) [Fig. IMaUb), Sec. 14.1.3. 2j . 

The more 9 initially decreases from tt the greater the smallest 2D Ra c (k y ) 
near k y = diverges from such J Ra c (fcy)|6»= 7r /| cos6*| [Fig. [TBT b). 9 > 1.5tt/2]. 
However, the neutralization of along-slot gravity by small G\ kz >o also makes 
the 3D analogue of the 2D branch with oj c {k y )\ ky ^Q x 1 to additionally connect 
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to that of the higher-o; c 2D branch with cu c (k y )\k y ->o — [Fig. I20l a) and (b)]. 
For k z > 0, closed contours of finite (Ra c , k y ) and (u c , k y ) thus form at these 
9. Growing with it — 6, such contours collide with the respective smaller 3D 
Ra c (k y ) and co c (k y ) around 9 ps lAir/2 (at the smaller n — 9 the smaller k z is, 
starting from the reorganization of the 2D branches). Via such a 3D collision, 
a structure with two connected limit points unfolds [Fig. 1207 c)]. Its hysteresis 
reconciles the small- and large-G small-k z branches of disparate 3D nature. 

For small k z , in particular, Ra c (k y , k z ) and u c (k y , k z ) are approximated by 
the Ra c (k y , 0) and u c (k y , 0) at the same 9. Such are the upper branches in 
Fig. I2"07 c). They fail to exist only for such small G (i.e., k y at a fixed k z ) as 
make the 3D effect of along-slot gravity relatively negligible. Largely triggering 
convection by the across-slot gravity alone, the small-G nearly-steady lower 
branches in Fig. 1207 c) fail to exist when the increase of G (i.e., of k y at a fixed 
k z ) transforms the convective pattern so that it ought to be driven by the 
along-slot gravity as well. Some quantitative details are given in Table [2J 

With k decreasing along the upper branch, the 2D effect of along-slot gravity 
grows with respect to that of across-slot gravity, according to the relative 
roles of the across-slot and along-slot density differences. Neutralization of 
such along-slot gravity in (3D) Eqs. ( 112"]) and (113"]) would thus require G to 
decrease with k at the upper limit point. With dG/dk > at this limit point, 
dk y /dk z > k y /k z , analogously to (j3"7|) following from (13T)]) . 

Along the (lower) branch largely specified by k alone, however, dG/dk > 
dG/dk > 0. Arising due to the across-slot gravity alone, the lower 

branch thus reaches its upper bound when G becomes so large that the latter 
inequality fails. At this bound then dG/dk < 0, as discussed in Sec. 4.1.2 of [32] 
for | tan#| and the exactly steady 2D analogue of such a branch at r\ — 1 and 
X = 0. For the lower limit point, therefore, dG/dk z < and dk y /dk z < k y /k z . 
This restriction leads to vanishing of the hysteresis region with increasing k z 
[Fig. 1201 c) and Table |2]. Such a region also vanishes as 9 approaches n/2. 
For k z (> 0) being fixed, in particular, Eq. fl38|) at both limit points implies 
dk y (k z ,6)/d6 > for the relevant 9 > vr/2 (Table [2]). For rj = 1 and % = 
(Fig. 10 in [32]), the exactly steady analogue of the current small-G branch 
also shrinks to k = as 9 reaches 7r/2. 

4.1.3.6 3D disturbances for rj = and x — !■ As discussed in Sec. 
l4.1.3.3l above. the relatively abrupt changes in 2D dRa c /dk y and du c /dk y for 
7] = and x = 1 are due to a switching between the primary roles of the two 
gravity components in the nature of the perturbation. They are thus specified 
by 9 alone. With k substituting for k y , the 3D effect of along-slot gravity 
decreases with respect to that of across-slot gravity, due to the emergence 
of G(< 1) in Eqs. ( fl2|) and (fl3|) . The 3D boundary at which the nature of 
perturbation changes, k y (k z ), thus deviates downwards from q y (k z ) = [k y (0) 2 — 



35 



k 2 ] 1 / 2 . This is seen from Fig. l2~TT a) and (b) and from Table [31 When the change 
in 2D dRa c /dk y is steep enough [Fig. l2~]T a)]. such a deviation leads to a short 
interval of k y with 3D most unstable disturbances separated from k y = 0. 

Conditions I and II are met for the steady marginal-stability curve at 9 = ir 
(9 = 0) [Figs. 171 and [Mf b)]. Here condition II holds due to the opposition 
provided by the effect of component conditions at the inverse (distinction) 
boundary to a steady cell rotation. For small ir — 9 (9) [Fig. l2TT a)]. 3D slow- 
propagating cells arising from the steady disturbances [localized at the dis- 
tinction (inverse) boundary] are thus most unstable near k y = 0. 

With Re[du(k,0)/dx]\ e=7T ^ oik ) x Re[du(k, 0)/dx)\ d=o ^ x u c (k, 0)| e=7r _ o(fc) = 
w c (fc,0)| e=o(fe) x 1 and Re[du(k,0)/dx]\ e=n / 2 x u c (k, 0) le=7r/2 x 1 as k -)> 

[Fig.LHb)], G|*-x> = 0(k)(^ 0) <=► [Rt o{k) {k)/G]\ k ^ = [R n i 2 o{k) (k)/G} lk ^ x 
[k/G]\k^Q > 0. Such a three-dimensionality then also exists so long as 9 ^ n/2 
[Fig. [211(b) and (c)], since the across-slot gravity is made dominant by a suf- 
ficiently small G at any 9 G (tt/2,iv) [9 G (0, 7r/2)]. This three-dimensionality 
could be viewed as coming from conditions I and II at G = as well. As 
k y — > for a fixed k z > 0, in particular, the Ra c (k y , k z ) \ cos9\ tends to the 2D 
Ra c {k z ) at 9 = n (9 = 0) with the u c (k y , k z ) ]ky ^ Q -> [Figs. [H(b), and EE]. 

Associated in Table [3]with ky(k z ), the relatively abrupt changes in dRa c /dk y 
and dojc/dky arising from k y = are thus also due to a transition between 
the perturbation suited mainly for the across-slot gravity and that for the 
along-slot gravity. Such a transition is however essentially of a 3D nature. 

With increasing k, the 2D effect of along-slot gravity decreases with respect 
to that of across-slot gravity. At ky(k z ) for a fixed 9 G (tv/2,tv) [9 G (0, 7r/2)], 
the latter gravity component still gives way to the former in driving the per- 
turbation. To neutralize the 2D across-slot gravity in Eqs. f[T2l and (fI3"]) . G 
must then grow with such k: dG/dk > for k = [k^{k z ) + k 2 ] 1 ' 2 . As (07} 
follows from flU}, therefore, dk G Jdk z > k^/k z > (Fig. |2T] and Table [3]). With 
i?e[<i£t(/c, O)/d:r]|0 =7r _ o (fc) x i2e[dtt(A;, 0)/c?x] |e =0 (^) x 1 and i?e[rf-u(A;, O)/rfx]|0 =7r / 2 x 
1 as k -> [fc 2 ->■ 0, fc„(A*) = ^(A; z ) = o(A; 2 )], [K-W^A?/] | fc ^ = 
[- R o+ (fe)(^)^/y|fe-^o - 1 implies (k 2 /k y )\ k ^ x 1 for fcj,(fc 2 ) = fcJ^A:*) [= o(fc*)]. 
Relations f[35l) then apply to k y (k z ). 

For large enough |0 — 7r/2| [Fig. l2TT a) and (b)], k y (0) is distinguishable within 
the considered range of k y . Since k 6 y {k z ) decreases with increasing k z , the 
above increase of the respective k y (k z ) with k z eventually leads to vanish- 
ing of the interval of k y where the instability is mainly due to the along- 
slot gravity. This takes place when k y (k z ) and k y (k z ) merge. When 9 is so 
close to it/2 that such k y {k z ) are not found [Fig. l2~lT c)]. the increases of 
Ra c (0,k z ) = Ra c (k z ) \e=o^/\ cos 9\ and of k y {k z ) itself with large enough k z 
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make the changes in dRa c /dk y and dco c /dk y at k y {k z ) indistinguishable. The 
slow-propagating cells localized near the distinction (inverse) boundary then 
smoothly transform into the traveling cells occupying the whole slot width. 

As discussed above, k 6 y {k z ) increases when \9 — ir/2\ decreases for a given 
k z . For k z (> 0) being fixed, however, Eq. (I3"g|) implies dky(k z ,9)/d9 > 
[dk£(k z ,9)/d6 < 0] for 9 e (n/2,ir) [9 e (0,vr/2)]. The ultimate vanishing of 
the relatively abrupt changes in dRa c / dk y and du c / dk y thus has to take place 
at the larger k z the closer 9 is to ir/2, as in Fig. [21] and Table [3j 

4-2 Linear steady instability for 9 = n 

4-2.1 Background 

Computed from Eqs. ([24|) and (I2"5"j) and boundary conditions ( TT4|) . (fTo']) . and 
f fl6|) . steady marginal-stability curves for 9 = ir and prescribed Ra are illus- 
trated in Figs. [221 [231 and [241 Their most unstable wave numbers, k c , and the 
respective Ra^(k c ) are provided in Table [H Although Fig. [24] and its data in 
Table [Hare also relevant to 9 = (as mentioned in the captions of Fig. [241 and 
Table S]), they are discussed in Sec. 14.2.2.31 below only in terms of 9 = ir. The 
shape of the presented curves is specified by the combined effects of across-slot 
diffusion and the first dissipation mechanism (Sec. 14. 1. 1.2)) . 

For Ra = 0, such diffusion is stabilizing. Since it is absent for x — 0, the most 
unstable wave number is zero, where the minimal effect of the first dissipation 
mechanism is achieved. Derivable via the long- wavelength expansion [4"5foT)] . 
the respective Ra s c (0) (Fig. [22] and Table [4J depends on the velocity conditions 
at the slot boundaries. 

For x = 1, the stabilizing effect of diffusion for Ra = renders the most 
unstable wavelength finite. Compared to the stress-free boundary condition, 
the enhancement of along-slot dissipation at a no-slip boundary reduces the 
across-slot portion of actively moving streamline particles with an along-slot 
density difference. Like diffusion, the effect of such disparity in along-slot dis- 
sipation is also augmented with the wavelength. Besides raising Ra s c (k) for all 
k, therefore, the no-slip effect increases k c as well (Figs. [23] and [24] in Table 
[4] for Ra = 0) [49.50J. The smaller across-slot portion of streamline particles 
driving convection then constitutes a larger part of the overall streamline. 

When introduced only at one boundary for Ra = 0, the no-slip effect is more 
pronounced in combination with the flux solute condition, where the solute 
perturbation isolines are orthogonal to the boundary [Figs. 123Tb) and [24T b). 
and Tabled]. This takes place because for the stress- free condition, the actively 
moving streamline particles with along-slot density differences exist infinites- 
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imally close to the boundary at which there is no diffusion. The (zero- .Ret) k c 
in Table @] (x = 1) for 7_ = 1 and 7+ = (7_ = and 7+ = 1) is thus closer 
to its value for j± = 1 (7± = 0) than to that for j± = (7± = 1) [32J50]. 



4.2.2 Ra>0 

4.2.2.1 i] = x = 0. When Ra increases from 0, the effect of diffusion 
becomes destabilizing. Being due to differential gradient diffusion at the dis- 
tinction boundary, the instability mechanism is then similar to that described 
in [28"f29"] . For 77 = x — 0, its effect is maximized at the infinite wavelength, 
where any stable fixed- value stratification is fully neutralized by diffusion. Also 
minimizing the role of the first dissipation mechanism, such wavelength thus 
remains most unstable for any Ra [49J. For k — > 0, the instability parameters 
in Fig. [22] (Table [4} are identical to those for rj = 1 and x — [28"|29|I32|I49|I51| . 

When k increases from 0, however, quantitative differences from the respective 
Ra s c {k) for 77 = 1 and x — arise, due to differential gradient diffusion taking 
place only at one boundary. In Fig. 122T b). such differential diffusion at the 
distinction boundary is also more effective in combination with the respective 
stress-free single boundary condition (7- = and 7+ = 1) than with the 
no-slip one (7- = 1 and 7+ = 0). As discussed above, the (nondiffusive) flux 
component perturbation is affected by the disparity between the no-slip and 
stress-free velocity conditions more than the (diffusive) fixed-value one. 



4.2.2.2 rj = x — !• When coupled with the distinction boundary, the 
single stress-free condition is more destabilizing than in its combination with 
the similarity boundary for 77 = x = 1 as well. In particular, this explains the 
relative location of the solid and dashed lines in Fig. [231( b). With growing Ra, 
in addition, the Ra s c (k c ) in Table H]for 7_ = and 7+ = 1 becomes even much 
closer to such a value for 7± = than to that for 7_ = 1 and 7+ = 0. Likewise, 
such a Ra s c (k c ) for 7_ = 1 and 7+ = becomes much closer to its respective 
value for 7± = 1 than to that for 7_ = and 7+ = 1. 

That the effect of diffusion turns destabilizing when Ra increases from is 
also manifested for = x = 1 in decreasing k c (Fig. [23] and Tabled]). Due 
to solute neutralization by its diffusion at the similarity boundary, however, 
k c cannot decrease to zero. It is the balance between the overall differential 
effect of both boundaries on component perturbations and the effect of the 
similarity boundary on solute perturbation that specifies the value of k c . 

As the overall effect of differential diffusion is augmented with growing Ra, 
however, the ratio of Ra s c (k) to Ra decreases for any k. This relatively enhances 
the effect of the similarity boundary on solute perturbation in the above bal- 
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ance specifying k c , at which Ra s c (k)/Ra is the smallest for a fixed Ra. When 
the latter effect turns dominant above certain Ra, k c begins to increase. As 
in Table HI such an oscillatory behavior of k c has to be most pronounced for 
7± = 0, where diffusion is not restricted by dissipation at the boundaries. It 
has to be least pronounced for j± = 1, where k c = 0.86 for Ra = 200000. 



4.2.2.3 r] = and x = 1. For 77 = and x — 1 (Fig- Ell), the increase of 
Ra from introduces differential gradient diffusion that favors both a steady 
perturbation at the distinction boundary and a standing-wave perturbation 
at the inverse boundary. Such a diffusion process at the latter boundary thus 
opposes a steady rotation of the convective cells. As suggested by Fig. [7J 
however, the steady linear instability for 7_ = 7+ is expected to precede the 
oscillatory one. (As indicated in the captions of Figs. 0, [181 and EH the results 
for 9 = and for 9 = 71 are relevant to each other.) Whereas both dissipation 
mechanisms accentuated in Sec. 14.1 .1.21 above damp the oscillatory (standing- 
wave) perturbation, only the first one of them damps the steady perturbation. 

Let the velocity conditions be either the same at both boundaries (7. = 7+) or 
such as the single stress-free condition is at the distinction boundary (7. = 
and 7+ = 1). Localization of steady convective cells near the distinction bound- 
ary is thus favored even when Ra is small. (For 8 = 0, such a localization at 
the inverse boundary is seen in Fig. [TBI) The destabilizing effect of the distinc- 
tion boundary then dominates the stabilizing effect of the inverse boundary. 
Thus k c only decreases with Ra increasing from for such 7^ (Table ED . 

When the single stress-free condition is at the inverse boundary (7_ = 1 and 
7+ = 0), however, the effect of differential gradient diffusion at the distinction 
boundary may not be dominant. At small Ra, in particular, the onset of 
steady convection is dominated by the disparity between the velocity boundary 
conditions. Indeed, the overall effect of diffusion on steady convection at Ra = 
1000 remains stabilizing: Ra s c (k c ) | Ra=1000 = 1949 > Ra+Ra s c (k c )^ Ra=0 = 1000+ 
817 i Table [Q). 

For 7_ = 1 and 7+ = 0, k c thus first grows with Ra increasing from (to 
at least 1000, as in Table H]). Only when the stable fixed- value stratification 
becomes large enough, the overall effect of differential gradient diffusion on 
steady convective cells becomes destabilizing. With the cells localized near 
the distinction boundary, in particular, k c decreases as Ra increases to 5000 
and above (Table ED- For all pairs of •y±, however, its decrease to is still 
prevented by the stabilizing effect of the inverse boundary. At the infinite 
wavelength, such effect becomes commensurate with the destabilizing role of 
the distinction boundary. 

Compared to the no-slip boundary condition, the stress-free condition en- 
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hances differential gradient diffusion at its boundary, as discussed above. Such 
an enhancement is also more relevant at the boundary where the convective 
cells are driven. In particular, the dashed lines are above the respective solid 
lines in Fig. 123T b) . For growing Ra, in addition, Ra s c (k c ) for 7_ = and 7 + = 1 
remains closer to such a value for 7± = than to that for 7_ = 1 and 7+ = 0. 
The Ra s c (k c ) for 7_ = 1 and 7+ = also remains closer to such a value for 
7± = 1 than to that for 7_ = and 7+ = 1 (Table H for Fig. EM]) . 

4-3 Finite- amplitude steady convection for 9 = 
4.3.1 77 = x= 1 

For the component conditions being different only at one boundary and 6 = 0, 
the most pronounced manifestation of the finite-amplitude steady instability 
mechanism reported for rj = 1 and \ = in (37] takes place when 7] = \ = 1. 
Its bifurcation diagram is illustrated in Fig. [25] When the relative role of stable 
solute stratification increases with \i, the linear steady instability is delayed for 
two reasons. One of them is merely the stabilizing solute contribution to the 
overall background stratification. The other reason is of a double-component 
nature. It is the opposition provided by differential gradient diffusion to a 
steady rotation of convective cells. Such finite-amplitude convection is however 
relatively little affected by the growing solute stratification, due to the nature 
of its mechanism. This nature is illustrated in Fig. [26] 

As the convection amplitude increases [Fig. 1267 a) .(b)], a growing number of 
solute isolines near the distinction boundary are found to be "outside" the 
flow domain. This effect is particularly pronounced in the regions of across-slot 
motion towards this boundary. The background solute scale is thus reduced 
in the convective state. Due to the fixed-value temperature conditions at both 
boundaries, however, the isotherms only increase their density in the direction 
of across-slot motion near either boundary. The background temperature scale 
thus remains intact even when convection is well-developed [Fig. 1257 b) — (e)]. 

Substantially reducing the across-slot solute scale, the finite-amplitude pertur- 
bation thus generates such unstable density gradients in the regions of across- 
slot motion [Fig. !2B7 a) — (e)] as exceed the linear single-component thresh- 
old [Ra c {n) ~ 1707 [52]]. This takes place well before onset of the double- 
component linear instability (Fig. 123]) . Horizontal density differences between 
the respective streamline points are then formed and give rise to (positive) 
convective feedback. Coming from finite-amplitude Rayleigh — Benard convec- 
tion, this feedback maintains the disparity between component gradients. 

Such finite-amplitude convection also arises when the overall background strat- 
ification is neutral or stable. For the present formulation, the linear steady 
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instability is then absent [Fig. l26T c) — (e)]. With limit point L in Fig. [251 mov- 
ing to higher Ra for growing /i, the flow amplitude and the disparity between 
component scales at the same Ra (= 60000) still decrease [Fig. l2"6T d).(e)]. 

Continuation of such a finite-amplitude steady flow for /z = lin6 l to#>0 
fails even when the step in 9 is close to zero. Such a failure of the contin- 
uation procedure is also experienced when an attempt to continue a super- 
critical finite-amplitude steady solution for // = 1 and 9 = n to 9 < n is 
made. The latter failure has to be due to the anticipated transformation of 
the steady finite-amplitude solutions into traveling-wave ones, as discussed for 
the respective small-amplitude flows say in Sec. 14.1.3.11 above. Purely finite- 
amplitude transformations of a traveling wave into a steady flow have been 
reported in binary-fluid convection [53] and in conventional double-diffusive 
convection [51]. The continuation failure at 9 = could thus also be due to a 
transformation of the steady finite-amplitude flow into a traveling-wave one. 



4-3.2 Comparison with 77 = 1 and x — and 77 = \ < 1 

Differential gradient diffusion at both boundaries delays the linear steady in- 
stability for 9 = more effectively than this process at the distinction bound- 
ary alone. In the finite-amplitude manifestation of steady convection, the en- 
hancement of differential diffusion on account of the second boundary also 
plays a stabilizing role. Solute stratification is however substantially reduced in 
such a convective steady state. Being partly due to differential gradient diffu- 
sion, its (stabilizing) effect on finite-amplitude steady convection is then much 
smaller than that on the respective small- amplitude one. When jj, < 1, the hys- 
teresis thus has to be more pronounced for 77 = 1 and x — than for f] = X = 1- 

In particular, the interval of hysteresis is larger for fi = 0.6 in Fig. 2 of [37] 
than for \i = 0.7 in the present Fig. [25J (see also Table |5J. Limit point L 
for 77 = x — 1 still precedes that for 77 = 1 and x — at the same fi. 
Compared to r] = x = 1, the stabilizing effect of differential gradient diffusion 
on finite-amplitude steady convection is enhanced for 77 = 1 and x = more 
than the efficiency of such finite-amplitude instability. [Indeed, the convection 
amplitude in Fig. l2"6T c) is slightly higher than that in Fig. 4(c) of [27J-] For 
this reason, limit point L for rj — % — 1 continues to precede that for 77 = 1 
and x = at least when fi e [1, 1.5], where the linear steady instability sets 
in at infinite Ra in both these cases. In particular, this is seen from Fig. 2 
(including the caption) in J27J and the present Fig. [25] for fj, — 1 and fi = 1.5. 

[Regarding the solute scale and phase specification, the present numerical 
formulation for 77 = x = 1 (Sec. 12. ip is not identical to that in [37]. Using 
the same steady solution phase and no other restrictions, however, it is most 
consistent with the formulation in [57] . The above relative location of limit 
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point L for n = \ = 1, with respect to that for n = 1 and x = 0, was 
also found to take place when the solute concentration values are fixed at the 
middle points of the slot boundaries in both cases.] 

The finite-amplitude mechanism discussed in Sec. 14.3.11 above is underlain by 
the flux condition decreasing in the perturbed state the across-slot scale of a 
component on which it is imposed at a boundary. The across-slot perturbation 
scales of both components then have to be substantially smaller for n = x = 
than for r) — x — 1- The resulting smaller disparity between such scales is 
thus generally less likely to trigger finite-amplitude steady convection before 
the onset of small-amplitude one. This explains why, for j± = 1, the hysteresis 
practically vanishes when f] = x approaches at least for /i < 0.9. 

When ?7 = x < 1 is sufficiently above 0, however, the hysteresis region persists. 
In particular, it is indisputably present at least for n = x ~ 0.41 (Tabled 7± = 
1). Such a finite-amplitude steady flow is illustrated in Fig. l2"oT f). As discussed 
in Sec. 14.2.21 above, differential gradient diffusion is also more effective at a 
stress-free slot boundary than at a no-slip one. Disparity between the small- 
and finite-amplitude manifestations of such a process is an important factor for 
the bifurcation subcriticality. The stress-free boundary conditions thus ought 
to have a quantitative effect on the hysteresis region. 

For fi < 1 (Table [5]), in particular, the subcriticality for n = x ~ 0.41 is 
more pronounced for 7± = than for 7± = 1. (The potentially interesting 
cases 7_ ^ 7+ are outside the scope of this work.) Trial computations with 
7 ± = at fi G [0.7, 0.9] also suggested that such a hysteresis then survives for 
rj = x = as well. However, the latter finding could not be ascertained because 
of a numerically transcritical manifestation of the hysteresis at n = x = 0. This 
ought to be due to an unfolding of the bifurcation diagram by the scale-fixing 
condition, to which the formulation for rj = x = is particularly sensitive. 

4-3.3 Environmental implications 

Let Ra in Fig. [25] be fixed on branch A2 above limit point L. This implies a 
value of /i where bifurcation point B is subcritical. For such a flow, let fi (and 
Ra s ) be varied from a negative value to the (positive) value just defined and 
above. Fig. [25]then suggests the existence of hysteresis in /i as well. One could 
thus also expect such a hysteresis to arise when Ra and Ra s are independent of 
each other. In particular, it ought to take place when Ra s is varied sufficiently 
for Ra being fixed on the higher- amplitude branch of the respective limit point. 

Viewed as a possible explanation of abrupt climate changes [22|23] , such a 
behavior has been reported in numerous model studies of the global ocean 
thermohaline circulation and climate [20|21f23p4^ . In these studies, the fresh- 
water flux into a North Atlantic ocean region of deep sinking (deep-water 
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formation) is a parameter controlling the salt stratification there. 

[Hysteresis was also found for salinity- and temperature-driven regimes of a 
box model considered in [TJ]. The boundary condition for a component is then 
specified by the component conductivity at the boundary Such hysteresis is 
viewed as a simplified illustration of multiplicity of the regimes of global ocean 
thermohaline circulation and climate [16 20.21.22 23 24]. The boxes in [15] can 
however represent the upper and lower ocean layers. The mechanisms of the 
salinity- and temperature-driven regimes [H] are then conceptually analogous 
to the hydrodynamic mechanisms of double-component convection in [2"8f2"5] 
(Sees. 14.2.2.11 and 14.2.2.21 above) and in [3?j (Sees. jZH and 12X21 above), 
respectively. 

In the framework of such analogy, in particular, the difference between box 
values of either component corresponds to the component gradient in the 
respective convective state and the rate of hydraulic flow to the convection 
amplitude. The external mechanical mixing is also represented by the (equal) 
component eddy diffusivities. (If one disregards the physical dissimilarity be- 
tween component boundary conductivities and their fluid diffusivities, such 
an analogy could also be respectively extended to the small- [TJ and finite- 
amplitude [2] steady instabilities in conventional double-diffusive convection.)] 

A basic constituent of the global ocean thermohaline circulation is comprised 
by the so-called Atlantic Meridional Overturning Circulation (AMOC). Op- 
eration of AMOC implies the existence of a North Atlantic region of deep 
sinking [TTf2 0,22,23j. Such a region also gives rise to penetration of horizon- 
tal sea-surface density differences into the ocean depth, without which AMOC 
could not be effectively driven by buoyancy [20,55,56,57]. In the context of [55] 
discussed in [2"0"|56|57j . therefore, the nature of an AMOC driven by buoyancy 
could only be such as suggests no source of the AMOC hysteresis other than 
the sinking region. The convection hysteresis, if any, in this region thus has to 
be a local factor imposing the solution structure and the parameters of hys- 
teresis on the entire AMOC (and the respective climate regimes {22123124] ). 

Allowing for transformation fllOp . however, a hysteresis between steady regimes 
in the convection region is certain if there is a hydrodynamic formulation 
where the temperature condition at a relevant boundary other than the sea- 
surface boundary is far enough from the flux type, as discussed above. Such a 
formulation could arise for the flow domain just above the bottom topography. 
Hysteresis in convection then ought to exist regardless of its parameters near 
V — X — for the stress-free boundaries if the eddy diffusivity (defined without 
convection) within the topographic bottom layer is not too small compared 
to that within the upper flow domain itself. This involves the issue of vertical 
variation of such an eddy diffusivity. With many uncertainties, its discussions 
[20,56,58j still seem to suggest that the lower boundary condition has to be 
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within the range of existence of a significant hysteresis. 



5 Summary and concluding remarks 

5. 1 General 

This work provides a comprehensive insight into the manifestation of double- 
component convection due to different boundary conditions in a diversely ori- 
ented infinite slot when the reflection symmetry between the slot conditions for 
a component is broken. In a class of the addressed problems, different compo- 
nent conditions at one slot boundary (the distinction boundary) are considered 
primarily with flux and with fixed-value conditions for both components at 
the other (the similarity boundary). Another class of the problems is such 
as the component conditions at the second boundary (the inverse boundary) 
differ from each other inversely to the distinction boundary. For elimination of 
other physical effects, the primary focus is on the compensating background 
gradients and equal component diffusivities. Mainly treating small-amplitude 
convection, this study also examines finite-amplitude steady instability. 

Being imposed only at one slot boundary, different component conditions still 
give rise to double-component convection for any slot orientation to the grav- 
ity. They also do so being inversely prescribed at the boundaries, although 
either component then has one fixed-value and one flux boundary condition. 
In either of these problem classes, however, the manifestation of small- and 
finite-amplitude convection is substantially disparate from that with the re- 
flectionally symmetric boundary conditions for a component. 

One aspect of such a disparity can be interpreted as coming from dissimilar 
2D small-amplitude convection patterns and the ranges of their formation. It 
is largely a consequence of only one traveling wave being relevant when the 
slot orientation differs from horizontal. Another cause of the disparity is a 
more frequent dominance of 3D small-amplitude disturbances, due to both a 
previously identified three-dimensionality mechanism and new ones. Directly 
underlying these factors, the broken symmetry also has major implications for 
finite-amplitude steady instability. 

5.2 Inviscid fluid 

When the similarity boundary is of the flux type (rj = x — 0), the behavior of 
an inviscid oscillatory marginal-stability curve in the horizontal slot for 6 = 
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is least dissimilar from that under the reflectionally symmetric component 
boundary conditions (77 = 1 and x = 0). In particular, such a linear stability 
boundary is characterized by the most unstable wave number being zero. For 
independently prescribed Ra s , Ra c and du c /dk at k = obtained from the 
long-wavelength expansion also differ from the symmetric case only by the 
denominator in the expression for the former parameter. Compared to 77 = 1 
and x — 0, the lack of temperature diffusion at the similarity boundary for 
77 = x = stabilizes most wavelengths. However, it also destabilizes such long 
wavelengths as effectively neutralize the temperature scale by diffusion at the 
distinction boundary alone. 

As at i] = 1 and x = 0, the zero wave number at rj = x — remains most 
unstable for the compensating background gradients up until the slot orienta- 
tion becomes nearly vertical. For 9 > 0, however, the reflection asymmetry of 
oscillatory perturbation leads to the expression relating 2D Ra c and du c /dk y 
at k y = being entirely dissimilar from that for rj = 1 and x — 0. In addition, 
zero limit of the Ra c {k y ) and ui c {k y ) as 9 — > 71/ 2 does not necessarily imply 
that the oscillatory marginal-stability boundary transforms into a steady one. 
For any 9 G [0,7r), no steady linear instability arises in viscous fluid at finite 
Ra c for either f] = X- The inviscid instability with Ra c = u c = for 9 > tt/2 
could not thus be underlain by viscous steady instability. This also applies to 
other such inviscid oscillatory-instability zero thresholds for 9 G (0, it). 

Despite the absence of dissipation, the most unstable wave number for the 
oscillatory instability at 9 = becomes finite both when 77 = x = 1 and when 
77 = and x — 1- F° r V — X = 1> such a change comes only from solute diffusion 
at the similarity boundary. For rj = and x — lj it is due to differential gradi- 
ent diffusion at the inverse boundary [or at the distinction boundary for 9 = n, 
allowing for (T2"Uj) ]. When the stable gradient is independently prescribed, ei- 
ther of such dissimilar diffusion processes still also prevents a manifestation 
of the oscillatory instability in an immediate vicinity of k = 0. 

For 77 = and x — 1 a t 9 = and 9 = n, however, the vanishing of oscillatory 
instability with decreasing k is also preceded by a short interval of k where a 
higher unstable gradient of one component destabilizes the lower stable strat- 
ification of the other. This is attributed to the growing relative disparity be- 
tween the oscillation frequencies at such different stable stratifications as these 
frequencies decrease with the increasing wavelength. Differential diffusion at 
the boundary favoring steady convection is thus enhanced relatively stronger 
by the lower stable stratification. At small enough k, the unstable gradient is 
then affected more by such stabilizing enhancement than by the effect of the 
stable stratification at the boundary favoring oscillatory convection. 

When 9 = in inviscid fluid, solute diffusion at the similarity boundary for 
77 = x — 1 allows to meet (the second of) two general conditions for three- 
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dimensionality of instability disturbances formulated in [32]. Due to these con- 
ditions, herein referred to as conditions I and II, such a three-dimensionality 
then arises in the vicinity of k y = for small 9 > 0. As 9 increases further 
to 7r/2, however, 2D Ra c (k y ) and u c (k y ) for the compensating background 
gradients decrease to zero when rj = \ = 1 as well. It is also the vicinity of 
k y = that is most sensitive to the direct effect of along-slot gravity, and 
thus experiences the fastest decrease of 2D Ra c (k y ) and co c (k y ). Unlike other 
3D instability manifestations resulting from conditions I and II, therefore, the 
present three-dimensionality is not retained above small 9 > 0. 

In contrast to the above problems for r\ = x, an interval of zero 2D Ra c (k y ) 
and u) c (k y ) arises for r\ = and x — 1 under the compensating background 
gradients already when 9 (tt — 9) is infinitesimal. [Invariance (|20|) also makes 
9 G (0, 7r/2) for r] = and x = 1 at /i = 1 equivalent to (tt — 9) G (0, 7r/2).] 
Expanding with growing 9 (tt — 9), this interval becomes infinite as 9 — > tt/2. 
As its lower limit also increases with 9 (tt — 9), however, the Ra c (k y ) and 
uj c (k y ) below such a limit are left finite for any 9 G (0, 7r/2] [9 G [tt/2,tt)). 
Neutralizing either component by diffusion, such relatively long wavelengths 
prohibit the 2D instability manifestation at infinitesimal Ra(k y ). 

The upper limit of the interval of zero 2D Ra c (k y ) and u> c (k y ) thus formed for 
r] = and x — 1 when 9 G (0,7r/2] [9 G [7r/2,7r)] is interpreted as resulting 
from a competition between two dissimilar patterns. Arising for the relatively 
longer wavelengths, one of them has steadily rotating cells driven by both the 
along-slot gravity and the effect of across-slot gravity at the inverse (distinc- 
tion) boundary. It is a traveling wave whose speed turns infinitesimal when 
Ra c (k y ) = u) c (k y ) = 0. The other pattern is featured by cells whose sense 
of rotation changes periodically in time. Characterized only by finite Ra c (k y ) 
and u c (k y ), it arises from the effect of across-slot gravity at the distinction 
(inverse) boundary for the relatively shorter wavelengths. 

Although conditions I and II are met in inviscid fluid at 9 = (9 = tt) for rj = 
and x = 1 as well, 3D instability then arises for 9 G (0,7r/2) [9 G (tt/2,tt)] 
largely from more pronounced effects. Two of these effects are due to the 2D 
interval with Ra c (k y ) = oo c (k y ) = 0. One of them is underlain by an invariance 
[Eq. (|3T|) ] of the linear stability equations at 9 = tt/2. With the 2D area of 
Ra c (k) = 0, this invariance leads to an interval of k y where 3D disturbances 
are most unstable for 9 = tt/2. Such an interval thus arises at small \tt/2 — 9\ 
as well. For any 9 G (0, tt/2) [9 G (tt/2,tt)], however, 3D disturbances are 
also dominant due to their unequal effect on the two gravity components. 
Relatively reducing the along-slot gravity by G = k y /k, this effect shifts the 
3D zero-threshold interval of k y below such (2D) interval specified by 9 alone. 

3D most unstable disturbances in inviscid fluid for rj = and x — 1 a ls° arise 
from another type of 3D areas with zero Ra c (k y , k z ) and u c (k y , k z ). Such a 3D 
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area has no 2D progenitor. Taking place at any 9 G (0, 7r/2) [9 G (vr/2, 7r)], such 
a perturbation three- dimensionality is still consistent with conditions I and II 
when they are considered at G — k y /k = 0. The nature of the zero-threshold 
area underlying such a 3D instability is interpretable in terms of the analogy 
between the effect of G in the 3D marginal-stability equations for k and the 
effect of |tan#| in the 2D equations for k y . With this analogy, differential 
behavior of the boundaries of such an area with respect to k z and 9 is consistent 
with the respective numerical data. It also matches the asymptotic behavior 
of the area boundaries for k z — >■ obtained by an independent approach. 



5.3 Viscous fluid 

For all the considered combinations of rj and Xi oscillatory linear instability 
also arises in viscous fluid. When 9 = 0, in particular, the marginal-stability 
curve for either 77 = x is characterized by an infinite increase of Ra c (k) with 
k decreasing to 0. Both such increases come from reduction of the instability 
feedback efficiency by along-slot dissipation. First identified for 77 = 1 and 
X = in [32J , this mechanism of efficiency reduction is herein referred to as the 
second dissipation mechanism. For 77 = x = 1, however, the long wavelengths 
are also so stabilized by solute diffusion at the similarity boundary. 

The viscous instability at 9 = for 77 = x — 1 still precedes that for 77 = x — 
so long as k is not too large: 8Ra c (k) > for such k. Compared to 77 = x — 0, 
diffusion of both components at the similarity boundary then leads to efficiency 
of the instability mechanism for 77 = x — 1 being higher. Such higher efficiency, 
however, generates steeper velocity gradients and thus more dissipation. Rel- 
atively increasing with k, this effect of dissipation eventually dominates the 
effect of diffusion at the similarity boundary. SRa c (k) thus changes its sign. 
With k increasing further, it also exhibits decaying oscillations. The oscilla- 
tory behavior results from an additional dissipation arising above a value of k 
where 5Ra c (k) = acting against the increase of \5Ra c (k)\. The additional dis- 
sipation is greater for such 77 = x as has the higher efficiency of the combined 
effects of diffusion and dissipation at the given k. 

The second dissipation mechanism is the reason why oscillatory linear insta- 
bility for rj = and x = 1 at # = and 9 = n is preceded by the steady one. 
For sufficiently long wavelengths, such oscillatory instability at 9 = is pre- 
ceded by those for rj = x as well. The stabilizing effect of the inverse boundary 
for 7] = and x = 1 is then more important than the cell localization near 
the distinction boundary. As the wave number increases, however, the role of 
this stabilizing effect diminishes. Compared to 77 = x, the cell localization for 
77 = and x = 1 a l so leads to a faster growth of the relative portion of fluid 
particles with along-slot density differences, and thereby to the faster enhance- 
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merit of efficiency of the instability mechanism. Above certain k, therefore, the 
oscillatory instability for r\ = and x = 1 at # = precedes those for rj = x- 

In viscous fluid, conditions I and II are met at 9 = for either rj = x- 111 
both cases, 3D disturbances are thus most unstable in a vicinity of k y = 
when 9 > is small. For either r] = x, however, 3D disturbances with small 
k y also remain the first to arise so long as 9 < tt/2. This takes place because a 
sufficiently small G = k y j k renders the effect of across-slot gravity dominant 
at any 9 G (0,7r/2). Compared to the inviscid fluid for r\ — x — 1, this is 
possible due to the respective viscous 2D effects of along-slot gravity being 
near k — sufficiently weak relatively to those of across-slot gravity. 

In viscous fluid, the broken symmetry between boundary conditions for a 
component is manifested most commonly at 9 — ir/2. Oscillatory instability 
is then characterized by a sequence of 2D counter-rotating convective cells 
propagating with a nonzero speed in the gravity direction. For any pair of 
considered i] and x-, t ne respective propagating pattern gives rise to such a 
distribution of the component perturbations as favors its convective motion 
and maintains the direction of its propagation. This direction is sustained by 
a nearly steady horizontal density difference generated between the traveling 
convective pattern and slot boundaries. Such a mechanism of traveling-wave 
instability is underlain by differential gradient diffusion. 

Although differential gradient diffusion for rj = and x = 1 ac t s inversely at 
the slot boundaries, its sidewall effects at 9 = it/ 2 match each other in the 
framework of the above features of the instability pattern. Such inversely sym- 
metric effects of differential diffusion at 9 = tt/2, however, result in horizontal- 
ly of across-slot convective motion, in contrast to i] = x- The corresponding 
reduction of the across-slot cell path then leads to a smaller relative con- 
vection amplitude of the instability. As a consequence, the marginal-stability 
frequencies are respectively lower for rj = and x — 1 than for r] = x- 

Horizontal density differences comprise the key quantitative feature of the 
vertical-slot instability mechanism. Their independence of orientation of the 
isolines of both components at the similarity boundary is thus the reason why 
the marginal-stability curves at 9 = n /2 are independent of 7] = x- The effect 
of along-slot dissipation on the efficiency of instability feedback at 9 = n/2 is 
also similar to that of the second dissipation mechanism at 9 — 0. Although 
the marginal-stability curve for r] = x = lat^ = 0is affected by component 
diffusion at the similarity boundary as well, this effect does not introduce a 
qualitative change. Its role also gradually vanishes as 9 = n/2 is approached. 
Continuous transformation of either 2D marginal-stability curve for f] = x 
thus takes place as 9 increases from to ir/2. 

The nature of vertical-slot instability to the single gravity-directed traveling 
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wave suggests a universal consequence of the broken symmetry in viscous fluid. 
In particular, such a traveling-wave character of the instability is maintained 
for any considered pair of r\ and x U P until 9 increases to tt, where the in- 
stability becomes steady. Any of these 2D transformation scenarios is also 
accompanied by a manifestation of 3D oscillatory disturbances. Such a mani- 
festation is underlain by different features of the respective steady instability 
at 9 = 7i. Embodied at 9 < tt in a small-frequency behavior of the respective 
3D oscillatory-instability boundary, these dissimilar features thus make the 
specific 3D scenario of transformation depend on the values of rj and x- 

The simplest 3D scenario of the above transformation takes place for rj = 
X = 1. It is specified by the nature of 3D disturbances being merely due to 
conditions I and II at 9 — tt. In this region of 9, however, such conditions imply 
that the instability is generally of an oscillatory nature with only 9 = tt being 
specifically characterized by oo c (k) = 0. The small value of u c > is thus a part 
of the perturbation introduced by the deviation of 9 from tt. An interval of k y 
where 3D disturbances are dominant is then retained by a properly decreased 
G = ky/k when 9 is close to tt/2 as well. This is permitted by the appropriate 
relative behavior of the respective 2D effects of the along-slot and across-slot 
gravity components near k — 0. 

For rj = x = 0, condition II is not met at 9 = tt. However, a vicinity of k y = 
is then substantially more unstable to 2D disturbances at 9 = tt than at 9 G 
(tt/2, tt). Mimicking at small G = k y / k the behavior of the 2D steady marginal- 
stability boundary at 9 = tt, such a 3D oscillatory-instability curve at 9 G 
(tt/2, tt) thus gives rise to 3D most unstable disturbances. With growing tt — 9, 
the large-G behavior of the small-A^ 3D linear-instability curve increasingly 
diverges from the small-G one. Isolated contours of solutions of the 3D linear- 
stability equations then also form out of secondary multiple solutions of such 
2D equations. The small- and large-G types of 3D behavior are thus reconciled 
via a hysteresis region between them. Such a region arises from collision of the 
isolated solution contours with the 3D boundary of primary instability upon 
the respective reorganization of the 2D marginal-stability solution structure. 

Such a strongly nonlinear behavior of solutions of the linear oscillatory-instability 
equations is an outcome of the broken symmetry and the resulting universality 
of oscillatory manifestation of the effect of different boundary conditions. For 
rj = x = and small tt — 9, such a universality transforms the solution of the 
linear steady-instability equations for i] = 1 and x — into a small-frequency 
solution of such equations for oscillatory perturbation. This furnishes condi- 
tions for the solution multiplicity of the latter 3D equations. Utilizing the 
physics of such multiple solutions, behavior of the resulting hysteresis region 
is interpreted based on the above analogy between the 3D effect of G = k y /k 
and the 2D effect of |tan#| in the respective linear stability equations. Such 
interpretation leads to the differential properties of this region with respect to 



49 



k z and 9 that are consistent with the numerical data. 

For r) = and x = 1, abrupt marginal-stability changes arise when 6 G (ir/2, it) 
[9 G (0,7r/2)]. One their class comes from such a 2D change. Its 2D nature 
is due to an irreconcilableness between the effects of along-slot and across- 
slot gravity. The failure of the former to accommodate the latter is related 
to a phase shift between the component and flow perturbations. Such a shift 
arises from the opposition generated by the across-slot gravity at the inverse 
(distinction) boundary to steady cell rotation. Above certain k y , the convection 
pattern thus switches from that largely driven by the along-slot gravity to a 
slower-propagating (like nearly-steady) one driven mainly by the across-slot 
gravity. The latter pattern is localized at the distinction (inverse) boundary. 
According to the sensitivity to the phase shift, the respective change in the 2D 
marginal-stability curve arises at the longer wavelength the larger \9 — n/2\ is. 

Due to the effect of G = k y /k(< 1) in the 3D linear-stability equations, the 
3D analogue of an above 2D abrupt change for rj = and x = 1 shifts to 
smaller k y as k z increases from 0. When such a 2D change is steep enough, as 
for \9 — tt/2\ > 37r/8, this shift leads to a small region isolated from k y = 
where 3D disturbances are most unstable. Another region of 3D most unstable 
disturbances arises then near k y = 0. Being due to conditions I and II at 9 = ir 
(9 = 0), it is maintained for any 9 G (n/2,n) [9 G (0, 7r/2)] by a sufficiently 
small G. This is allowed by virtue of an appropriate weakness of the 2D effect 
of along-slot gravity with respect to that of across-slot gravity near k = 0. 

3D marginally unstable disturbances arising when 9 G (7r/2,7r) [9 G (0, 7r/2)] 
for rj = and x=\ from k y — exhibit another class of abrupt changes. Being 
essentially of a 3D nature, these changes however also imply switching between 
a pattern driven by the across-slot gravity and that driven by the along- 
slot gravity. The former is characterized by slow-propagating convection cells 
localized at the distinction (inverse) boundary while the latter by traveling 
convection cells occupying the whole slot width. Differential behavior of such 
a transition with respect to k z and 9 is described with the above analogy 
between the 3D effect of G and the 2D effect of |tan#| in the respective 
marginal-stability equations. 

Steady linear instability at 9 = n arises for any pair of r/ and \ due to the 
mechanism of differential gradient diffusion identified in [28 29J. When the 
slot boundaries have dissimilar velocity conditions and rj — Xi such a mech- 
anism triggers the instability for the stress-free condition at the distinction 
boundary before that for the no-slip one. This takes place because the stress- 
free condition accommodates the formation of along-slot density differences 
infinitesimally close to its boundary. Such an effect could have manifestations 
beyond the instability for which it was exposed. Although differential gra- 
dient diffusion for 77 = and x — 1 a t 9 = n (6 = 0) acts (inversely) at 
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both boundaries, the steady perturbation is driven at the distinction (inverse) 
boundary. The steady instability for the single stress-free velocity condition 
at this boundary thus still precedes that for the single no-slip one. 

For rj = Xi t ne distinction boundary alone could also give rise to finite- 
amplitude steady instability at 9 = 0. Compared to r] — 1 and x — [37], such 
an instability is less subcritical for \i < 1. Its onset at r\ = x = 1 still precedes 
that for rj = 1 and x — 0, for the second boundary in |37j is also a source 
of stabilizing differential diffusion. When both such linear steady instabilities 
set in at infinite Ra with /i > 1, the hysteresis for rj = x — 1 is thus even 
more prominent. Still well-pronounced at r\ = x ~ 0.41, it practically vanishes 
for the no-slip boundaries at rj = x — 0. The convective scales of both com- 
ponents are then reduced too much for their disparity to give rise to such a 
subcriticality. The hysteresis at rj = x < 0-41 could also be more prominent 
for stress-free slot boundaries. The possibility that abrupt climate change be 
largely controlled by such a finite-amplitude instability is thus discussed. 



5.4 Main conclusions 

Some consequences of the broken symmetry are associated with 2D small- 
amplitude convection. At 9 = 0, in particular, component diffusion at the 
similarity boundary for r] = x = 1 is found to enhance efficiency of the viscous 
standing- wave instability compared to rj = x = 0. The key outcome in viscous 
fluid is however the universally oscillatory manifestation of small- amplitude 
convection. Such a universality stems from the preference given by the broken 
symmetry to one of the two counter-traveling waves over the other. At 77 = 
and x — 1) it also involves an abrupt change in the 2D marginal-stability 
curve for an inclined slot, when dissimilar oscillatory patterns are switched. In 
inviscid fluid, such a change comes with an interval of 2D wave numbers where 
Ra c = bj c = 0. Also arising at 9 = n/2, such an interval is then a separate 
cause for three-dimensionality of most unstable disturbances near 9 = 71/ 2. 

The broken symmetry also generally enhances the role of 3D small-amplitude 
disturbances. Some of them arise from conditions I and II. Such a mechanism is 
also relevant when the nature of instability changes from steady to oscillatory 
with the parameter value at which conditions I and II hold being perturbed. 
Besides the second dissipation mechanism, the 3D effect of condition II arises 
from solute diffusion at the similarity boundary for ^ = x = 1, both at 
9 = and at 9 = 7T, and from differential gradient diffusion at the inverse 
(distinction) boundary for rj = and x = 1 at ^ = 7r (^ = 0). For 9 = 0, 
the former case is also relevant to inviscid fluid. Unlike the other 3D effects 
coming from conditions I and II, such 3D inviscid disturbances at 77 = x = 1 
are dominant only for small 9 > 0. This results from the 2D along-slot gravity 
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near k — being then relatively strong with respect to the across-slot gravity. 

3D disturbances also arise from or involve mechanisms other than that spec- 
ified by conditions I and II. For r\ = and % = 1, their dominance when 
9 G (0, 7r/2) [9 G (n/2,7i)] is associated with abrupt changes in the 3D 
marginal-stability curves both in inviscid and in viscous fluid. In inviscid fluid, 
in particular, such abrupt changes arise with Ra c = u c = 0. In viscous fluid 
for r) = x — 0, such a three-dimensionality stems from the disparity between 
the 2D behavior of Ra c near k y = at 9 = it and that at 9 G (it/2,it). This 
disparity also leads to multiplicity and isolated existence of as well as hystere- 
sis between solutions of the 3D linear oscillatory-instability equations. The 
3D behavior of the inviscid and viscous abrupt marginal-stability changes for 
i] = and x = 1 as well as of the hysteresis region for rj = \ — follows from 
the analogy between the 3D effect of G — k y /k and the 2D effect of | tan#|. 

The broken symmetry also has other important implications. In the context of 
linear steady instability at 9 = ir, it exposes a higher effectiveness of differential 
gradient diffusion at a stress-free fluid boundary than at a no-slip one. It also 
has a two-fold effect on the mechanism of finite-amplitude steady instability 
at 9 — 0. Eliminating for rj = \ one boundary at which gradient disparity 
could form for r\ — 1 and x = 0, such a symmetry breaking also reduces the 
stabilizing effect of differential gradient diffusion on account of this boundary. 
The convection hysteresis then effectively arises for r\ — x — 1- Although 
such a hysteresis could practically vanish when the convective scales of both 
components are reduced at r\ — x — 0, it remains significant when rj = x is 
closer to than to 1. Interpretation of abrupt climate change as being largely 
underlain by such a finite-amplitude convective instability is thus discussed. 
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o-^ = 5s - 5t 
(l-r|)t x +r|t=(r|-2) R a/2 (l-x)s x +%s=(%-2)Ra^2 

! / if 

t = Ra/2 s x = - Ra s 

Fig. 1. Geometry of the problem. 5p = p — po is the variation of the (dimensionless) 
density, p, due to the variations Ss and 6t of solute concentration s and temperature 
t with respect to their reference values, at which the density is po; a = gd 3 /nv. 
Pr = u/k = 6.7, Le = kt/ks = 1; kt and «s(= «) are the component diffusivities. 
The fluid is of the Boussinesq type. 
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Fig. 2. 9 = 0. Inviscid fluid. Curves of the marginal linear stability to 2D oscillatory 
disturbances for independently prescribed Ra s , Ra c (k) and ui c (k); Le = 1. The solid 
lines: Ra s = 50000, the dashed lines: Ra s = 20000, the dotted lines: Ra s = 10000, 
the dash-dot lines: Ra s = 5000. (a) rj = x = 0; (b) rj = x = 1; ( c ) V = 0, X = 1 [hi 
view of transformation ([20]) . these data apply to 9 = ir as well if Ra c \-> Ra s c and 
Ra s !->■ i?a]: in particular, Ra c (k) for iia s = 5000 exceeds that for Ra s = 10000 when 
k < 0.56, Ra c (k) for Ra s = 10000 exceeds that for Ra s = 20000 when k < 0.37, 
and Ra c (k) for i?a s = 20000 exceeds that for Ra s = 50000 when k < 0.25. 
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Fig. 3. Inviscid fluid; rj = x = 0- Curves of the marginal linear stability to 3D 
(k z > 0) oscillatory disturbances for different k z , Ra c (k y ) and u c (k y ); fj, = 1, Le = 1. 

(a) k z = 0, the solid lines: 9 = 0, the dashed lines: 6 = 7r/8, the dotted lines: = vr/4; 

(b) k z = 0, the solid lines: 9 = 37r/8, the dashed lines: 6 = 0.97r/2, the doted lines: 
9 = 0.99vr/2; (c) 6 = 0.05vr/2, the solid lines: k z = 0, the dashed lines: k z = 0.5, the 
dotted lines: k z = 1, the dash-dot lines: k z = 2. 
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Fig. 4. Inviscid fluid; rj = x = 1- Curves of the marginal linear stability to 3D 
(k z > 0) oscillatory disturbances for different k z , Ra c (k y ) and u c (k y ); \i = 1, Le = 1. 
For (a) — (e), the solid lines: k z = 0, the dashed lines: k z = 0.5, the dotted lines: 
k z = 1, the dash-dot lines: A; z = 2. (a) = 0; (b) = 0.05vr/2; (c) 9 = O.lvr/2; (d) 
9 = 0.15tt/2; (e) 9 = tt/8; (f) k z = 0, the solid lines: = tt/4, the dashed lines: 
9 = 37r/8, the dotted lines: 9 = 0.9tt/2, the dash-dot lines: 9 = 0.99vr/2. 



60 




Fig. 4. Inviscid fluid; 77 = % = 1. Curves of the marginal linear stability to 3D 
{k z > 0) oscillatory disturbances for different k z , Ra c {k y ) and co c (k y ); /j, = 1, Le = 1. 
For (a) — (e), the solid lines: k z = 0, the dashed lines: k z = 0.5, the dotted lines: 
k z = 1, the dash-dot lines: k z = 2. (a) = 0; (b) 9 = 0.05vr/2; (c) 6 = O.br/2; (d) 
(9 = 0.15tt/2; (e) = tt/8; (f) A; z = 0, the solid lines: 9 = tt/4, the dashed lines: 
9 = 3tt/8, the dotted lines: 9 = 0.9tt/2, the dash-dot lines: 9 = 0.99vr/2. 
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Fig. 5. Inviscid fluid; r\ = 0, % = 1. Curves of the marginal linear stability to 
2D (k z = 0) oscillatory disturbances, Ra c (k y ) and uj c (k y ); (j, = 1, Le = 1. The 
nearly vertical curves of Ra c (k y ) are expected to go to Ra c (k y ) = 0. The solid 
lines: 9 = (1 =f 1)tt/2, the dashed lines: = (1 T 0.99)tt/2, the dotted lines: 
9 = (1 =F 0.95)vr/2, the dash-dot lines: 9 = (1 ^ 0.9)vr/2. The two values of 6» are 
relevant in light of transformation (I20j) . 
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Fig. 6. Inviscid fluid; r] = 0, % = 1. Curves of the marginal linear stability to 
3D (k z > 0) oscillatory disturbances for different k z , Ra c (k y ) and uj c (k y ); fi = 1, 
Le = 1. The nearly vertical curves of Ra c (k y ) are expected to go to Ra c (k y ) = 0. (a) 
= (1=f0.85)vt/2; (b) 9 = (1t0.75)tt/2; (c) 9 = (1^0.5)vr/2; (d) 9 = (1 T0.25)vr/2; 

(e) 9 = (1 =F 0.1)7r/2; (f) 9 = tt/2. The two values of 9 in (a) — (e) are relevant in 
light of transformation ([20]) . The solid lines: k z = 0; the dashed lines: (a) k z = 0.15, 
(b) k z = 0.25, (c) k z = 0.5, (d) k z = 1.3, (e) k z = 1.9, (f) k z = 1; the dotted lines: 
(a) k g = 0.4, (b) k z = 0.6, (c) k z = 1.5, (d) k z = 2.2, (e) k z = 3, (f) k z = 1.5; the 
dash-dot lines: (a) k z = 0.5, (b) k z = 0.77, (c) k z = 2, (d) k z = 2.8, (e) k z = 4.2, 

(f) k z = 2. Additional data are provided in Table HJ 

63 




02460246 
k k 

y y 

Fig. 6. Inviscid fluid; r] = 0, % = 1. Curves of the marginal linear stability to 
3D (k z > 0) oscillatory disturbances for different k z , Ra c (k y ) and uj c (k y ); fi = 1, 
Le = 1. The nearly vertical curves of Ra c (k y ) are expected to go to Ra c (k y ) = 0. (a) 
= (1=f0.85)vt/2; (b) 9 = (1t0.75)tt/2; (c) 9 = (1^0.5)vr/2; (d) 9 = (1 T0.25)vr/2; 

(e) 9 = (1 =F 0.1)7r/2; (f) 9 = tt/2. The two values of 9 in (a) — (e) are relevant in 
light of transformation ([20]) . The solid lines: k z = 0; the dashed lines: (a) k z = 0.15, 
(b) k z = 0.25, (c) k z = 0.5, (d) k z = 1.3, (e) k z = 1.9, (f) k z = 1; the dotted lines: 
(a) k g = 0.4, (b) k z = 0.6, (c) k z = 1.5, (d) k z = 2.2, (e) k z = 3, (f) k z = 1.5; the 
dash-dot lines: (a) k z = 0.5, (b) k z = 0.77, (c) k z = 2, (d) k z = 2.8, (e) k z = 4.2, 

(f) k z = 2. Additional data are provided in Table HJ 
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Fig. 7.9 = 0. Viscous fluid; rj = 0, x = lj A 4 = 1; -f r = 6.7, Le = 1. In light of 
transformation (120p . these results apply to # = ir as well. Curves of the marginal 
linear stability to 2D steady (the solid lines) and oscillatory (the dashed lines) 
disturbances, Ra c (k) and u; c (k), for stress- free (s) and no-slip (n) slot boundary 
conditions. Notations s and n refer to the respective closest curves below them. 
(The analogue of the standing-wave oscillatory instabilities described in this figure 
is not reported for the inclined slot with viscous fluid considered herein below. It is 
then preceded by two types of 2D traveling-wave instability, one of which is related 
to the steady instabilities described in this figure.) 
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Fig. 8. 8 = 0. Viscous fluid and no-slip slot boundaries; fx = 1, Pr = 6.7, Le = 1. (a) 
Curves of the marginal linear stability to 2D oscillatory disturbances, Ra c {k) and 
uj c (k); the solid lines: 77 = % = 1, the dashed lines: 77 = x = 0. (b) The dotted lines: 
5Ra c (k) = Ra c (k) lr)=x=0 - Ra c (k) \ v=x=1 and 6u c (k) = uj c (k) lv=x=0 - u c (k)\ v=x 
the solid lines designate the respective zero values. 
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Fig. 9.8 = 0. Viscous fluid and no-slip slot boundaries; r/ = x = 0. Perturbation 
temporal behavior throughout a half of the oscillation period t p ~ 185t (St = 0.05) 
just beyond the onset of 2D oscillatory instability. It was obtained from the numer- 
ical simulation of evolution of the linearized Eqs. ([I]) — ([3]) in response to the initial 
disturbance proportional to the background state after initial time Tj ~ 340000 has 
passed; A = 2, /j, = 1, Ra = 87340, Pr = 6.7, Le = 1. With this Tj, all perturba- 
tion modes other than the unstable mode (r p ~ 18<5r) are practically negligible, tp: 
perturbation streamlines; s: isolines of solute concentration perturbation; t: pertur- 
bation isotherms. The actual relative values of the streamfunction perturbation are 
equal to 10 -3 times the respective values in the figure. The solid and dashed stream- 
lines designate the clockwise and counterclockwise rotation and are equally spaced 
within the positive and negative streamfunction intervals, respectively. The solid 
and dash-dot isolines of the component perturbations are equally spaced within the 
positive and negative component perturbation intervals, respectively, (a) r = Ti+5r; 
(b) T = n + 35t; (c) T = n + 55t; (d) T = n + 75t; (e) T = n + 8<5r; (f ) r = n + 95t. 
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Fig. 10. 8 = 0. Viscous fluid and no-slip slot boundaries; 77 = x = 1. Perturbation 
temporal behavior throughout a half of the oscillation period r p ~ 185t (St = 0.05) 
just beyond the onset of 2D oscillatory instability. It was obtained from the numer- 
ical simulation of evolution of the linearized Eqs. ([I]) — ([3]) in response to the initial 
disturbance proportional to the background state after initial time Tj ~ 120000 has 
passed; A = 2, \i = 1, Ra = 85040, Pr = 6.7, Le = 1. With this n, all perturba- 
tion modes other than the unstable mode (t p ~ 18<5r) are practically negligible. ij;: 
perturbation streamlines; s: isolines of solute concentration perturbation; t: pertur- 
bation isotherms. The actual relative values of the streamfunction perturbation are 
equal to 10~ 3 times the respective values in the figure. The solid and dashed stream- 
lines designate the clockwise and counterclockwise rotation and are equally spaced 
within the positive and negative streamfunction intervals, respectively. The solid 
and dash-dot isolines of the component perturbations are equally spaced within the 
positive and negative component perturbation intervals, respectively, (a) r = Tj+(5r; 
(b) r = n + 35t; (c) T = n + 55t; (d) T = Ti + 7<5r; (e) r = n + 8<5t; (f ) r = n + 95t. 
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Fig. 11. Viscous fluid and no-slip slot boundaries; f] = x = 0. Curves of the marginal 
linear stability to 3D (k z > 0) oscillatory disturbances for different k z , Ra c (k y ) and 
uj c (k y ); jj, = 1, Pr = 6.7, Le = 1. The solid lines: fc 2 = 0, the dashed lines: A: 2 = 0.5, 
the dotted lines: k z = 1, the dash-dot lines: k z = 2. (a) 9 = 0.25vr/2; (b) 9 = 0.5vr/2; 
(c) 6 = 0.75tt/2. 
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Fig. 12. Viscous fluid and no-slip slot boundaries; rj = x = 1- Curves of the marginal 
linear stability to 3D (k z > 0) oscillatory disturbances for different k z , Ra c (k y ) and 
uj c (k y ); jj, = 1, Pr = 6.7, Le = 1. The solid lines: fc 2 = 0, the dashed lines: k z = 0.5, 
the dotted lines: k z = 1, the dash-dot lines: k z = 2. (a) 9 = 0.25vr/2; (b) 9 = 0.5vr/2; 
(c) 6 = 0.75tt/2. 
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Fig. 13. Viscous fluid and no-slip slot boundaries. Curves of the marginal linear 
stability to 3D (k z > 0) oscillatory disturbances for different k z , Ra c (k y ) and u: c (k y ); 
H = l, Pr = 6.7, Le = 1. The solid lines: k z = 0, the dashed lines: k z = 0.5, the 
dotted lines: k z = 1, the dash-dot lines: k z = 2. (a) 9 = ir/2, rj = % = and 
K) = x = 1; (b) V = X = 0; (c) 7] = x = 1- hi (b) and (c), # = #7r/2, where 
the values of $ are given in these figures next to their respective closest curves the 
numbers do not intersect and for -d = 2, the instability is steady: uj c (k y ) = 0. For 
the lowest u c {k y ) > in (b) and for i9 = 2 in (c), the data near k y = (where such 
stable data were numerically difficult to obtain) are not presented. For such data, 
Ra c (k y ), k — > oo are assumed. 
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Fig. 14. Viscous fluid and no-slip slot boundaries; 77 = 0, % = 1. Curves of the 
marginal linear stability to 3D {k z > 0) oscillatory disturbances for different k z , 
Ra c (k y ) and u c (ky); fj, = 1, Pr = 6.7, Le = 1. The solid lines: k z = 0, the dashed 
lines: k z = 0.5, the dotted lines: k z = 1, the dash-dot lines: k z = 2. (a) 6 = tt/2; (b) 
6 = i?7r/2 and 8 = (2 — $)-7r/2, where the values of 1? are given in the figure next 
to their respective closest curves the numbers do not intersect and for # = 2, the 
instability is steady: oj c (k y ) = 0. The two values of 9 in (b) are applicable in light 
of transformation fj20f) . For 1? = 2 in (b), the data near k y = (where such stable 
data were numerically difficult to obtain) are not presented. 
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Fig. 15. = tt/2. Viscous fluid and no-slip slot boundaries; r/ = \ = 0. Perturba- 
tion temporal behavior throughout nearly a half of the oscillation period r p ~ 46<5r 
(5t = 0.05) just beyond the onset of 2D oscillatory instability. It was obtained from 
the numerical simulation of evolution of the linearized Eqs. ([!]) — ([3]) in response 
to the initial disturbance proportional to the background state after initial time 
n ~ 12000 has passed; A = 2, fi = 1, Ra = 11292, Pr = 6.7, Le = 1. With this n, 
all perturbation modes other than the unstable mode (r p ~ 46<5r) are practically 
negligible, ip: perturbation streamlines; s: isolines of solute concentration pertur- 
bation; t: perturbation isotherms. The actual relative values of the streamfunction 
perturbation are equal to 10~ 3 times the respective values in the figure. The solid 
and dashed streamlines designate the clockwise and counterclockwise rotation and 
are equally spaced within the positive and negative streamfunction intervals, respec- 
tively. The solid and dash-dot isolines of the component perturbations are equally 
spaced within the positive and negative component perturbation intervals, respec- 
tively. Small asymmetries between the positive and negative phases are due to the 
distortion introduced by the solute scale-fixing conditions. The instability onset Ra 
in the data underlying this figure is different from that for Fig. [16] only because the 
scale-fixing conditions make the numerical formulation used for these data slightly 
different from that used for the data in Fig. [T6J (a) r = Tj + 5t; (b) r = Tj + 6<5r; 
(c) t = n + 115t; (d) T = n + 1G5t; (e) r = n + 21<5r. 
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Fig. 16. = tt/2. Viscous fluid and no-slip slot boundaries; 77 = x = 1. Perturba- 
tion temporal behavior throughout nearly a half of the oscillation period r p ~ 46<5r 
(St = 0.05) just beyond the onset of 2D oscillatory instability. It was obtained from 
the numerical simulation of evolution of the linearized Eqs. ([1]) — ([3]) in response 
to the initial disturbance proportional to the background state after initial time 
n « 1200 has passed; A = 2, y, = 1, Ra = 11226, Pr = 6.7, Le = 1. With this n, 
all perturbation modes other than the unstable mode (r p ~ 46<5r) are practically 
negligible, ip: perturbation streamlines; s: isolines of solute concentration pertur- 
bation; t: perturbation isotherms. The actual relative values of the streamfunction 
perturbation are equal to 10~ 3 times the respective values in the figure. The solid 
and dashed streamlines designate the clockwise and counterclockwise rotation and 
are equally spaced within the positive and negative streamfunction intervals, respec- 
tively. The solid and dash-dot isolines of the component perturbations are equally 
spaced within the positive and negative component perturbation intervals, respec- 
tively. The instability onset Ra in the data underlying this figure is different from 
that for Fig. [15] only because a slightly different numerical formulation is used, (a) 
r = n + 5t; (b) r = 73 + 6<5t; (c) r = 73 + 115r; (d) r = t\ + 165r; (e) r = T{ + 21<5r. 
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Fig. 17. 9 = 7r/2. Viscous fluid and no-slip slot boundaries; rj = 0, % = 1. Perturba- 
tion temporal behavior throughout nearly a half of the oscillation period r p ~ 515r 
(<5t = 0.05) just beyond the onset of 2D oscillatory instability. It was obtained from 
the numerical simulation of evolution of the linearized Eqs. ([T]) — ([3]) in response 
to the initial disturbance proportional to the background state after initial time 
Tj ~ 1100 has passed; A = 2, \i = 1, Ra = 5586, Pr = 6.7, Le = 1. With this Tj, 
all perturbation modes other than the unstable mode (r p ~ 51<5r) are practically 
negligible, ip: perturbation streamlines; s: isolines of solute concentration pertur- 
bation; t: perturbation isotherms. The actual relative values of the streamfunction 
perturbation are equal to 10~ 3 times the respective values in the figure. The solid 
and dashed streamlines designate the clockwise and counterclockwise rotation and 
are equally spaced within the positive and negative streamfunction intervals, respec- 
tively. The solid and dash-dot isolines of the component perturbations are equally 
spaced within the positive and negative component perturbation intervals, respec- 
tively, (a) t = n+ 25t; (b) r = n + 75t; (c) t = r, + 13<5r; (d) r = n + 18<5r; (e) 
T = n + 235t. 
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Fig. 18. 6 = 0. Viscous fluid and no-slip slot boundaries; rj = 0, x = 1- P n light of 
transformation (|2Up . these results also apply to 9 = it if Ra c i— >• i?aj!, i?a s i-> i?a, 
and thus /U i— > l//i.] 2D singular eigenvectors corresponding to the wavelength A = 2 
at the onset of small-amplitude steady convection; Le = 1. ip: perturbation stream- 
lines; s: isolines of solute concentration perturbation; t: perturbation isotherms. The 
variables are nondimensionalized as in Eqs. (PQ) — ([3]). The actual relative values of 
the streamfunction perturbation are equal to 10~ 3 times the respective values in the 
figure. The solid and dashed streamlines designate the clockwise and counterclock- 
wise rotation and are equally spaced within the positive and negative streamfunction 
intervals, respectively. The solid and dash-dot isolines of the component perturba- 
tions are equally spaced within the positive and negative component perturbation 
intervals, respectively, (a) fi = 1, Ra c = 9787; (b) /i = 1.3, Ra c = 25478; (c) = 1.5, 
Ra c = 46908. 
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Fig. 19. Viscous fluid and no-slip slot boundaries; rj = x = 1- Curves of the marginal 
linear stability to 3D (k z > 0) oscillatory disturbances for different k z , Ra c (k y ) and 
Lo c (k y ); jJL = 1, Pr = 6.7, Le = 1. The solid lines: fc 2 = 0, the dashed lines: A; 2 = 0.5, 
the dotted lines: k z = 1, the dash-dot lines: k z = 2. (a) 6 = 1.757r/2; (b) = 1.57r/2; 
(c) = 1.25vr/2. 
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Fig. 20. Viscous fluid and no-slip slot boundaries; ry = x = 0- Curves of the marginal 
linear stability to 3D (k z > 0) oscillatory disturbances for different k z , Ra c (k y ) and 
u c (k y ); fj, = 1, Pr = 6.7, Le = 1. The solid lines: k z = 0, the dashed lines: k z = 0.2 
in (a) and (b) and k z = 0.5 in (c), the dotted lines: k z = 0.5 in (a) and (b) and 
k z = 1 in (c), the dash-dot lines: k z = 1 in (a) and (b) and k z = 2 in (c). (a) 
6 = 1.75-7T/2; (b) 6> = 1.5vr/2; (c) 6 = 1.25vr/2. The lowest-uj c (k y ) data near A; y = 
(where such stable data were numerically difficult to obtain) are not presented. 
Additional quantitative details are given in Table [2j 
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Fig. 21. Viscous fluid and no-slip slot boundaries; rj = 0, x = 1- Curves of the 
marginal linear stability to 3D {k z > 0) oscillatory disturbances for different k z , 
Ra c (k y ) and u c (k y ); fi = 1, Pr = 6.7, Le = 1. (a) 6 = (1 ± 0.75)vr/2; (b) 
(9 = (1 ± 0.5)vr/2; (c) 9 = (1 ± 0.25)vr/2. The two values of are relevant in light 
of transformation ()20p . The solid lines: = 0; the dashed lines: (a) k z = 0.5, (b) 
k z = 1, (c) k z = 2; the dotted lines: (a) k z = 1, (b) k z = 2, (c) k z = 3; the dash-dot 
lines: (a) k z = 2, (b) k z = 3, (c) k z = 4. Additional quantitative details are given in 
Table [5J 
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Fig. 22. 9 = ir. Viscous fluid; r/ = \ = 0- Curves of the marginal linear stability 
to steady disturbances, Ra s c {k), for different Ra; Le = 1. The numbers multiplied 
by 10 3 give the values of Ra for which the respective closest lines below them were 
computed, (a) the solid lines: j± = 0, the dashed lines: j± = 1; (b) the solid lines: 
7_ = and 7+ = 1, the dashed lines: 7_ = 1 and 7+ = 0. The most unstable wave 
numbers are described in Tabled! 
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Fig. 23. 9 = tt. Viscous fluid; 77 = % = 1. Curves of the marginal linear stability 
to steady disturbances, Ra s c {k), for different i?a; Le = 1. The numbers multiplied 
by 10 3 give the values of Ra for which the respective closest lines below them were 
computed, (a) the solid lines: j± = 0, the dashed lines: j± = 1; (b) the solid lines: 
7_ = and 7+ = 1, the dashed lines: 7_ = 1 and 7+ = 0. The most unstable wave 
numbers are described in Tabled! 
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Fig. 24. 9 = 7T. Viscous fluid; 77 = 0, x = 1. [These data also apply to = if 
transformation ([20]) along with Ra (->• i?a s , -Ra;! h-> i?a c , and 7± h-> 7^ are allowed 
for.] Curves of the marginal linear stability to steady disturbances, Ra s c {k), for 
different Ra; Le = 1. The numbers multiplied by 10 3 give the values of Ra for 
which the respective closest lines below them were computed, (a) the solid lines: 
7± = 0, the dashed lines: j± = 1; (b) the solid lines: 7_ = and 7+ = 1, the dashed 
lines: 7_ = 1 and 7+ = 0. The most unstable wave numbers are described in Table HI 

82 



\ * o 1^=0.7 |i=0.5 |i=0.3 |i=0.2 

A2 _ ^ V ^ =^ 

l _ . ' ;^ ^ 

g'-r - - - ->. (1=0.2 

4->X^<- "> 

B | - 

- ^= - 7 

B M 

(1=1.3 
(1=1.5 

tt 10000 20000 

Ra 




Fig. 25. = 0. Viscous fluid and no-slip slot boundaries; r) = x = 1- Schematic struc- 
tures of the 2D steady flows with minimal along-slot period A = 2 for /i € [0, 1.5]; 
Pr = 6.7, Le = 1. <& is an abstract measure of the steady flows that distinguishes 
between different solutions, specifies the location of the singularities (limit points 
and symmetry-breaking bifurcations), and represents the flows arising from a sym- 
metry-breaking bifurcation as a single branch. The background states are depicted 
by the horizontal lines with arrows (for fi = 1.5, this is the coordinate axis). The 
solid lines stand for the solutions being stable to the disturbances associated with 
the eigenvalues that give rise to the steady instability of the conduction state. The 
dashed lines represent the flows being unstable to either steady or both steady and 
oscillatory disturbances. The dotted lines stand for the solutions being unstable to 
oscillatory disturbances alone. Secondary bifurcations, if any, are not shown. B is 
the symmetry-breaking bifurcation standing for the steady linear stability boundary 
for wave number k = n (A = 2). Its criticality changes at \x just below 0.3. L is 
the limit point. Al and A2 are the lower- and higher-amplitude branches associated 
with the limit point, respectively. H is a Hopf bifurcation. For fj, = 1, it arises just 
below Ra « 85040 (Fig. HQ). 
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Fig. 26. 9 = 0. Viscous fluid and no-slip slot boundaries; Pr = 6.7, Le = 1, A = 2. 
2D convective steady flows, ip: streamlines; s: isolines of solute concentration; t: 
isotherms. The solid and dashed streamlines represent the clockwise and counter- 
clockwise rotation and are equally spaced within the positive and negative stream- 
function intervals, respectively. The actual values of s and t are equal to 10 3 times 
the respective values in the figure, (a) — (e): rj = \ = 1- ( a ) A 4 = 0.7, Ra = 9107, 
branch Al; (b) /j, = 0.7, Ra = 8000, branch A2; (c) /j, = 1, Ra = 24000, branch 
A2; (d) n = 1.3, Ra = 60000, branch A2; (e) /i = 1.5, Ra = 60000, branch A2; (f) 
r) = X ~ 0.41, /x = 1, ifa = 300000, branch A2. 
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Table 1 

Inviscid fluid; 7/ = 0, x = 1; // = 1, Le = 1. k 2l (k z ), k 2u (k z ), ky l (k z ), and ky U (k z ) are 
approximate values of wave number k y at which Ra c (k y ) and oj c (k y ) are expected 
to decrease to zero for different 9 and k z . k 2l (k z ) and k 2u {k z ) originate at k z = 0, 

whereas ky l (k z ) and ky U (k z ) arise only at > 0. For k z > 0, q l y (k z ) = [A^(0) - 

k 2 ] 1 ^ 2 and qy(k z ) = [k 2u (0) 2 — k 2 ] 1 ^ 2 are provided for comparison with k 2l (k z ) and 
k 2u (k z ), respectively. The two values of # are relevant due to transformation ([20]) . 



5) 


k 


k 31 
y 


l,3u 


k 21 

K y 


k 2u 

K y 


a 1 


u 
Qy 




n n 


Mot ton no 


M ot toi inn 


U.OuUU 










n 1 

U. x 


01 OQ 


031 9 


31 89 


8844 


3354 

U . OOOt: 


8898 




n 1 ^ 

U. XO 


094Q 


0840 

U.UO^X-U 


9880 

U.ZiUUU 


D 8704 

U.O ( U^i 


^1 89 

U.O-LUZ 


8897 


h x(l SKW/9 

I -L — p U. OO 1 /t / Zi 


n 9 


04^3 

U .UtiOO 


Nlot Ton nri 


Not ton no 
IsUL 1U 1111 U. 


8496 




8798 




u.o 


110^ 

U. 11UO 


M ot toi inn 


M ot TOI 1 n 
INUt 1UUX1U. 


7847 

U. / O^i 1 




84*^8 

U.O^iOU 




n 4 


9397 


Nlot Ton nri 


Not ton no 

1UUX1U. 


8899 




801 1 

0. ou ± ± 




^ 

U.O 


Mot Tfinnn 


M ot toi l n r\ 


M ot TOI 1 n 


M ot toi 1 n 








u.u 


1NUL 1UUIXU. 


1MUI 1UUIXU. 


u.oc/oy 


1 4^89 


u.oyoy 


1 4^89 




o 1 

u. 1 


o 0083 

u.uuuo 


01 8Q 

U.UXUa 


^78Q 

U.O / Oa 


1 44Q^ 


^874 
u.oo / *± 


1 4^98 




o 9 
u.z 


n 09^4 


07^3 
U.U / 00 


^908 
U.OZUU 


1 A303 


^81 ^ 
U.0D10 


1 4494 




n 9^ 

U . ZO 


040 1 

U.u L rU X 


1 399 


483Q 


1 41 ^4 


^40Q 


1 4348 




u.o 


0^8^ 


9741 

U.Zi / 'XX 


3991 


1 ^Qfifi 
i.oyuu 


^14Q 

U.O L'-kiJ 


1 49^0 






1 07Q 

U. 1U I u 


M ot toi i n H 


M ot TOI 1 n 


1 ^48^ 




1 4001 




u . o 


1 780 

O. X ( ou 


\Iot ton no 


\fnt ton no 


1 9781 




1 3677 

X . OU 1 I 




8 
u.u 


9784 


M ot toi inn 


M ot toi 1 n 
lift 1UUX1U. 


1 1 741 




1 3988 




u. / o 


8097 

U.UUZi / 


M ot toi inn 


M ot toi 1 n 


848^ 

U.04:00 




1 9489 




n 78 

U. 1 u 


8Q1 


Nlot ton no 


\fnt ton no 


7800 

U. ( UuU 




1 9491 




n 77 


Not mnnn 


Not toi i n n 


Nlot Ton no 


Nlot ton no 








0.0 


Not found 


Not found 


1.3006 


2.7783 


1.3006 


2.7783 




0.5 


0.0665 


0.2019 


1.1006 


2.7039 


1.2007 


2.7329 




0.6 


0.0971 


0.3266 


0.9756 


2.6697 


1.1539 


2.7127 


(1t0.5)tt/2 


0.7 


0.1344 


Not found 


Not found 


2.6284 




2.6887 




1.0 


0.2954 


Not found 


Not found 


2.4523 




2.5921 




1.5 


0.9178 


Not found 


Not found 


1.7937 




2.3386 




1.55 


1.0848 


Not found 


Not found 


1.6222 




2.3057 




1.6 


Not found 


Not found 


Not found 


Not found 







85 



Table 1 
(Continued) 
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Table 2 

Viscous fluid and no-slip slot boundaries; rj = x = 0; fi = 1, Pr = 6.7, Le = 1. ky{k z ) 
and k y (k z ) are approximate locations of the limit points of the upper and lower 
branches, respectively, of Ra c (k y ,k z ) and uj c (k y ,k z ) for different k z > and such 8 
as the limit points are connected to intervals of the primary-instability boundary 
[Fig. [201( c)] . Sign ps implies that the accuracy with which the number was obtained 
is smaller than that of the other numbers, due to very small values of uj c near such 
a limit point. For = (1 + 0.35)-7r/2, the limit points vanish just below k z = 1. 
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Table 3 

Viscous fluid and no-slip slot boundaries; 77 = 0, x = 1> A 4 = 1> P r = 6.7, Le = 1. 
ky{k z ) and k e y {k z ) are approximate values of wave number A;^ for different k z at 
a given where relatively abrupt changes in dRa c /dk y are distinguishable. (Such 
respective changes in dui c /dk y , in particular in Fig. I21L could be most pronounced at 
slightly shifted values of k y .) ky(k z ) originate at k z = 0, whereas ky(k z ) arise only 

at k z > 0. For k z > 0, q y{k z ) = [ky(0) 2 — k^] 1 ^ 2 are provided for comparison with 
ky(k z ). For 8 = (1 ±0.5)7r/2, ky{k z ) stands for the values of k y where the respective 
relatively abrupt change in doj c /dk y takes place at a given k z , since such a change 
in dRa c /dk y is then little distinguishable. The two values of 6 are relevant due to 
transformation (|20l) . 



9 


k z 




k e 
y 


4 


(1 ±0.75)7r/2 


0.0 


Not found 


2.3 


2.3 


0.5 


0.107 


2.2 


2.24 


1.0 


0.5 


1.7 


2.07 


2.0 


Not found 


Not found 


1.14 


(1±0.5)tt/2 


0.0 


Not found 


4.15 


4.15 


0.5 


0.04 


4.10 


4.12 


1.0 


0.175 


3.90 


4.03 


2.0 


0.8 


3.20 


3.64 


3.0 


Not found 


Not found 


2.87 


(1± 0.25)7r/2 


0.0 


Not found 


Not found 




0.5 


0.0175 


Not found 




1.0 


0.07 


Not found 




2.0 


0.3 


Not found 




3.0 


0.7 


Not found 




4.0 


1.6 


Not found 




5.0 


Not found 


Not found 
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Table 4 

9 = it. Viscous fluid; Le = 1. Approximate values of the most unstable wave number, 
k c , and of the respective Ra s c {k c ) for the linear steady instability at the values of 
Ra and j± used for Figs. [22], [231 and [Ml For 77 = and x = 1, these data also apply 
to = if transformation (|20p along with Ra 1— > Ra s , Ra s c 1— >• .Ra c , and 7± h4 7^ 
are allowed for. 








1000 


5000 


10000 


20000 


50000 


100000 


7? = x = 
[Fig. [22]^)] 


7± = 


k c 























Ra s c (k c ) 


120 


120 


120 


120 


120 


120 


120 


7± = 1 


k c 























Ra s c (k c ) 


720 


720 


720 


720 


720 


720 


720 


7? = x = 
[Fig. Elb)] 


7- = 0, 
7+ = 1 


k c 























Ra s c (k c ) 


320 


320 


320 


320 


320 


320 


320 


7 _ = 1, 
7+ = 


k c 























Ra s c (k c ) 


320 


320 


320 


320 


320 


320 


320 


V = x = 1 

[Fig. (Ha)] 


7± = 


kc 


1.76 


1.36 


0.97 


0.85 


0.77 


0.79 


0.85 


Ra s c (k c ) 


385 


882 


2539 


4392 


7825 


17007 


30186 


7± = 1 


k c 


2.55 


2.21 


1.54 


1.28 


1.09 


0.91 


0.85 


Ra s c (k c ) 


1296 


1980 


4284 


6744 


11167 


22973 


40590 


*7 = X = 1 

[Fig. EH1»] 


7- = 0, 
7+ = 1 


k c 


2.09 


1.71 


1.22 


1.04 


0.91 


0.82 


0.83 


Ra s c (k c ) 


669 


1183 


2870 


4704 


8048 


17024 


30241 


7_ = 1, 
7+ = 


kc 


2.21 


1.84 


1.27 


1.07 


0.93 


0.83 


0.84 


Ra s c (k c ) 


817 


1482 


3717 


6159 


10630 


22706 


40574 


^ = 0, x = 1 

[Fig. Ela)] 


7± = 


kc 


1.76 


1.75 


1.69 


1.59 


1.44 


1.22 


1.07 


Ra s c {k c ) 


385 


1289 


4003 


6437 


10246 


19327 


32200 


7± = 1 


kc 


2.55 


2.51 


2.20 


1.96 


1.76 


1.51 


1.32 


Ra s c (k c ) 


1296 


2255 


5745 


9445 


15500 


29334 


47784 


V = 0) X = 1 

[Fig. cap,)] 


7 _ = 0, 
7+ = 1 


kc 


2.09 


1.94 


1.68 


1.59 


1.49 


1.31 


1.16 


Ra s c (k c ) 


669 


1431 


4008 


6584 


10718 


20322 


33526 


7 _ = 1, 
7+ = 


kc 


2.21 


2.35 


2.20 


1.96 


1.72 


1.43 


1.23 


Ra s c {k c ) 


817 


1949 


5733 


9372 


15004 


27897 


45606 
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Table 5 

6 = 0. Viscous fluid; A = 2, Pr = 6.7, Le = 1. Obtained from the numerical steady 
solutions of Eqs. (pD) — ([3]), approximate values of Ra characterizing limit point L 
and symmetry-breaking bifurcation point B (Fig. [25]) for r/ = \ = 1 & t 7± = 1 an d 
such singularities for some = x < 1 at 7± = 1 and ■j± = in the range of \i < 1. 



A* 


0.3 


0.4 


0.5 


0.6 


0.7 


1 


»7 = X = 1 
7± = 1 


L 


2726 


3261 


3872 


4581 


5417 


9071 


B 


2733 


3419 


4547 


6694 


11834 


oo 


r] = x ~ 0.64 
7± = 1 


L 


Not found 


Not found 


3682 


4855 


6563 


25107 


B 


2299 


2844 


3716 


5300 


8800 


oo 


7/ = x ~ 0.41 
7± = 1 


L 


Not found 


Not found 


Not found 


4983 


7713 


269057 


B 


2194 


2708 


3524 


4992 


8184 


oo 


jy = X » 0.41 

7± = 


L 


Not found 


Not found 


1616 -7- 1617 


2274 


3343 


96701 


B 


980 


1223 


1616 -r 1617 


2342 


3938 


oo 
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